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FOREWORD 


This report presents, in graphical form, the results of studies of 
the eharacteristics of ballistic interplanetary trajectories to Mercury 
(launch dates: 1967-1968) and Jupiter (launch dates: 1968-1973). 
Also included are a description of the physical model, the development 
of the equations of the model, and a discussion of the properties of 
the trajectories. 
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ABSTRAQ 

The general characteristics of ballistic interplanetary trajectories are 
discussed, and detailed equations are developed for the analytical 
model. Extensive data are presented in graphical form for trajectories 
to Mercury (launch dates: 1967-1968) and Jupiter (launch dates: 
1968-1973). These graphs include: (1) curves of vis viva geocentric 
energy vs launch date for minimum-energy trajectories and (2) curves 
of 19-21 different trajectory parameters vs launch date for various 
vis viva geocentric energies. The trajectories were computed on the 
IBM 7090 digital computer by numerical evaluation of the analytical 
model, after which specific parameters of interest were automatically 
plotted, carefidly checked, and analyzed. Procedures are outlined for 
use of these data hy the trajectory engineer in the design and analysis 
of interplanetary trajectories. 


I. INTRODUCTION. ANALYTICAL MODEL FOR INTERPLANETARY TRAJECTORIES 


The analytical model used in the generation of Mer- 
cury and Jupiter trajectory parameters consists of three 
distinct phases of two-body motion: (1) an escape hyper- 
bola near the launch planet, (2) elliptical' motion under 
the attraction of the Sun, and (3) terminal hyperbolic 
motion near the target planet. 

A. He//ocenfr/c Alofion 

Solution of the heliocentric elliptic motion is obtained 
first under the following assumptions: 

( 1 ) The launch and target planets move in orbits about 
the Sun as given in the national ephemerides. 
Their velocity components are obtained by using 
two-body conic formulas, mean orbital elements, 
and their tabular positions, as listed in the ephem- 
erides. 

(2) The launch and target planets are massless. Thus, 
the only force acting on the probe is that of the 
Sun. 

( 3 ) The position of the probe at launch into the helio- 
centric orbit is the center of the massless launch 
planet. Its position at arrival on the heliocentric 
orbit is the center of the massless target planet. 


Thus, for solution of the heliocentric phase of motion, 
the attractions of the launch and target planets are tem- 
porarily disregarded. The primary result to be obtained 
from the solution of the heliocentric-transfer problem is 
determination of the hyperbolic-excess velocity vector 
relative to the launch planet. 


1. Determination of Planar Orientation 

Since the launch and arrival positions of the probe 
are assumed to be the centers of the launch and target 
planets, they can immediately be determined, given the 
launch and arrival^ times, by consulting the ephemeris. 
Further, the orientation of the heliocentric transfer plane 
can immediately be found. Let R/, be the Sun-launch- 
planet position vector at launch time T/„ and let Rp be 
the Sun-target-planet position vector at arrival time Ta 
(Fig. 1-1). Planar orientation is then found from the 
unit normal W to the plane, as follows: 


W 


^ Rr. X Rp 

Rl Rp sin ^ 


( 1 - 1 ) 


where the angle ^ is defined below. The inclination i to 


^Hyperbolic heliocentric motion is not considered herein. 


'Or, for convenience, the launch date and flight time can be specified. 
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the ecliptic plane^ can be found by 
cos i — W • K' 


(1-2) 



Fig. 1-1. Heliocentric-transfer geometry 


and the path angle T is found from 


^ = [V (1 - e^H2a - W ] ^ ^ ^ 

In Eqs. 1-6 and 1-7, GMs is the universal gravitational 
constant times the mass of the Sun ( =2.959122083 X 10"^ 
au^/day-); a and e are the semimajor axis and eccentricity 
of the transfer ellipse, respectively; and v is the true 
anomaly of the probe, given by 


cos V 


a{l - e^) -R 
eR 


(1-8) 


Now, there are two unknowns in Eqs. 1-5 to 1-8 which 
prevent their immediate evaluation: the semimajor axis 
a and the ( ccen tricity e. The determination of these 
quantities is the main problem. Battin (Ref. 1) has shown 
that the eccentricity is actually a function of the semi- 
major axis, rhus, it is first necessary to determine a. 
The semiinajor axis is related to the time of Hight Tr by 
Lambert’s Theorem, which states: The transfer time be- 
ticcc 7 i any two points on an ellipse is a fund km of the 
sum of the lUstances of each point from the focus, the 
distance hetueen the points, and the semimajor axis of 
the ellipse. Functionally, the theorem is stated as 


where K' is a unit vector pointing in the direction of 
the ecliptic north pole. 

2. In-Flane Relations 

The heliocentric central angle ^ (Fig. 1-1) is also 
readily determined by utilizing the positions of the launch 
and target planets. This angle may be obtained from 

“'■^ = 1511 ^ 

sin ^ = sgn [(Rl X Rp) • K'] (1 “ cos^ (1-4) 


Tr = TriRc + Rp.C,a) (1-9) 

where the distance C between the launch planet at 
launch time and the target planet at arrival time, as 
shown in Fig. 1-2, is obtained from 

C= |Rp-Rl| (1-10) 

Since the time of flight Tf and the launch and arrival 
positions R[, and Rp are known, only the semimajor axis 
remains to be found by iterative solution of Eq. 1-9. 
After the semimajor axis a is obtained, the heliocentric 
velocities of the probe at launch and arrival times V/. 


The velocity vector V of the spacecraft anywhere 
along its path may be obtained from 

v = ^[(W X R)cosr+ RsinT] (1-5) 
R 

Here, R is the heliocentric position vector, R = |R|, 
and V is the heliocentric speed, obtained from 

( 1 - 6 ) 



•In this Report, the interest is only in transfers which have the same 

rotational motion about the Sun as do the planets: thus, O^t P*9* 1—2, In-plane— transfer geometry 
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and Vp may be evaluated from Eq. 1-5 under the condi- 
tions R = Rl and R = Rp, The path angles F/,, Fp and 
true anomalies Vl, Vp at launch and arrival times^ may 
also be evaluated from Eqs. 1-7 and 1-8 under the same 
conditions. 

Finally, the desired end result, the hyperbolic-excess 
velocity \hr. relative to the launch planet, may be found 
(Fig. 1^3) by 

= V^- V, (1-11) 

where Vi is the velocity of the launch planet at launch 
time. 



Fig. 1-3. Determination of hyperbolic-excess 
velocity vector 


B. Launcfi-P/anef Escape Hyperbola 

The key result from the solution of heliocentric trans- 
fer is the hyperbolic-excess velocity vector VhL at launch. 
The reason for the importance of this vector is that it 
tells the direction in which the probe must be traveling 
relative to the launch planet when on the point of leaving 
the planet’s gravitational influence. There are an infinite 
number of escape trajectories ( all hyperbolas ) which can 
have the same hyperbolic-excess velocity vector. How- 
ever, only a portion of these are practical for use when 
related to existing launch sites and boost-vehicle con- 
straints. For example, it would be ridiculously costly in 
payload — and impractical — to shoot a vehicle straight up. 
Criteria for selection of a family of feasible escape tra- 
jectories are given below. 

1, Assumptions 

The solution of the escape phase of motion is obtained 
under the following assumptions: (1) The probe is acted 
on only by the gravitational force of the launch planet, 
and ( 2 ) the oblateness effects of the launch planet are 
neglected. 

^The details of quadrant choice for these angles are found in Ref. 3. 


The direction of the asymptote of the escape hyperbola 
is found by normalizing the hyperbolic-excess vector 
VftL- The injection energy C 3 of the escape hyperbola® is 
found by squaring the hyperbolic-excess speed, or 

C3 = V?^ (1-12) 

Thus, in contrast to the heliocentric problem, the launch 
planet is now “massy,” whereas the influence of the Sun 
is neglected. However, the hyperbolic-excess velocity 
vectors found by solving the heliocentric problem are 
used as a starting point to solve the escape problem. 

2. Size and Shape of Escape Hyperbola 

As previously stated, only some of the infinite number 
of escape trajectories are practical. Two of the practical 
aspects of a set of trajectories are the sizes and shapes 
of the hyperbolas. 

Size is basically determined by the energy C 3 , which, 
in turn, is a function of boost- vehicle capability. For 
boost vehicles in use at this writing (or shortly to be 
available), values of energy less than or equal to 25 
km^/sec^ are considered reasonable. The larger the value 
of energy that the booster is required to deliver, the 
smaller the payload and launch period over which the 
vehicle may be fired. 

The shape of the hyperbola is determined by its eccen- 
tricity, which is a function of both the energy and the 
perifocal distance, according to 


where Rp is the perifocal distance and GM is the uni- 
versal gravitational constant times the mass of the launch 
planet. From Eq. 1-13, it can be seen that, for a fixed 
perifocal distance, the eccentricity increases linearly with 
the energy. The value of perifocal distance is not arbi- 
trary, but depends strongly on the boost-vehicle trajec- 
tory. It has been shown (Ref. 2) that, in the great 
majority of cases, it is necessary and desirable to use a 
circular parking orbit as part of the preinjection phase of 
the escape trajectory. It is, further, an interesting fact 
that the altitude of the parking orbit determines the 
perifocal distance. If h is the parking-orbit altitude and 
Ro is the launch planet’s radius, then, to an extremely 
close degree of approximation, 

Rp = Ro + h (1-14) 

“Cais actually twice the total energy per unit mass; i.e., the uts ciua 
integral. 
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or the perifocal distance is equal to the launch-planet- 
centered radius of the parking orbit. In Ref. 2, it also 
has been shown that the lowest possible parking orbit 
(80-100 nm) allows greatest payload capability. Thus, 
using 100 nm for the parking-orbit altitude, a practical 
value of perifocal distance is 6560 km. The perifocal 
distance will vary only slightly about this value for other 
parking-orbit altitudes, or even for direct-ascent-type 
preinjection trajectories. Therefore, both the size and 
shape are essentially determined by the energy alone, 
which is found from Eq. 1-12. 

Given the size and shape of the escape hyperbola, its 
planar orientation must be determined. This can be done 
by considering two vectors: (1) the direction of the 
hyperbolic-excess vector, denoted by a unit vector S, and 
(2) a unit vector directed from the center of the 
launch planet to the launch site. The vehicle’s flight plane 
will essentially be determined by these two vectors, as 
shown in Fig. 1-4. A unit normal W to the launch-planet- 
centered flight plane is determined by 

X S 

'^“TRixsr 

with the constraint that the Z component of W is always 
positive. 

Since R[ is a function of time, according to the rotation 
rate of the launch planet, the planar orientation must 
continually change. In effect, this says that the launch 
azimuth is a continuous function of launch time. 

A detailed description of the geometrical aspects of 
the launch-planet ascent trajectory is not given here, but 
may be found in Ref. 2. 

C. Ca/cu/afion of Differential Corrections 

The calculation of differential corrections for inter- 
planetary trajectories may be accomplished in several 
ways and depends on the choice of independent and de- 
pendent variables. In this Report, a numerical differenc- 
ing scheme is used. Basically, the independent variables — 


the injection energy C3, declination and right ascen- 
sion 0s of the outgoing asymptote S of the escape hyper- 
bola — are varied, one at a time, to produce variations in 
the dependent variables — the components of the impact 
parameter B and the time-of-flight Tf. 

The impact parameter B is defined as a vector origi- 
nating at the center of the target planet and directed 
perpendicular to the incoming asymptote of the target- 
centered approach hyperbola (Fig. 1-5). The impact pa- 
rameter B is resolved into two components which lie in a 
plane normal to the incoming asymptote S. The orienta- 
tions of the reference axes in this plane are arbitrary, but 
one is usually selected to lie in a fixed plane. Thus, define 
a unit vector T, lying in the ecliptic plane, according to 

^-|sincT ^ ^ ^ 

where K' is a unit vector normal to the ecliptic plane. 
The remaining axis is then given by a unit vector R, 
defined by 

R = SXT (1-17) 

Figure 1-6 illustrates the orientation of the R, S, T 
target coordinates. 

The impact parameter B lies in the R-T plane and has 
miss coinponents B * T and B • R. The condition B*T~ 
B • R — 0 denotes vertical impact on the target. Thus, 
B • T, B • R, and Tr are the three target dependent vari- 
ables. If represents a set of generalized independent 
variables, such as injection position and velocity or other 
convenient variables, then the partial derivatives 
ciB • T/SQi, r'B • R/0Qi, dTjr/cQi are first-order differen- 
tial corrections or error coefficients relating miss at the 
target and flight- time errors to the independent variables. 

A convenient set of independent variables for inter- 
planetary trajectories is the vis viva injection energy C3, 
the declination and the right ascension 0,s- of the 
asymptote of the escape hyperbola. These variables essen- 
tially describe the launch hyperbolic-excess velocity vec- 
tor VftL. 
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Fig. 1^. Vehicle flight plane 
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II. DETAILED EQUATIONS FOR TRAJECTORY COMPUTATIONS 


In Section I, a summary was given of the physical 
model used to generate the parameters of Mercury and 
Jupiter trajectories. The purpose of this Section is to pre- 
sent, in full detail, the equations of the let Propulsion 
Laboratory’s Heliocentric Conic Trajectory Program^', in- 
cluding the equations of Lambert’s Theorem^ which were 
developed for this investigation. The equations were 
coded for the Laboratory’s IBM 7090 digital computer. 

This program has proved very useful over the past few 
years in studying and analyzing interplanetary trajectories. 
Some of the results of these studies are presented graphi- 
cally in Sections III to X of this Report. The program 
has been devised for calculation of trajectories from any 
planet in the solar system to any other planet. In addition, 
by adding the appropriate orbital elements, trajectories 
to comets or to any body of known motion may be 
computed. 

The Heliocentric Conic Trajectory Program is divided 
into four major sections: 

( 1 ) The heliocentric phase, in which basic computa- 
tions are made for a unique conic trajectory, of 
given flight time, that passes through the center 
of two massless planets rotating around the Sun, 
as defined in the national ephemerides. 


(4) Computation of partial derivatives which relate 
miss and time-of -flight errors at the target to vari- 
ations in the hyperbolic-excess velocity vector of 
the escape hyperbola near the launch planet. 

A. Heliocentric Phase 

1. Determination of Position at Launch and Arrival 
For any given time, the positions of the planets can be 
obtained from the ephemerides, referenced to a given 
coordinate system and epoch. For the present purposes, 
the heliocentric ephemerides of the planets, referenced 
to the mean equator and equinox of 1950.0, were selected. 

From these ephemerides, one may find the position vec- 
tor R' of the launch planet at launch time T/ , the posi- 
tion vector R' of the target planet at arrival time T,, and 
the position vector R' , of the launch planet at arrival 
time. For practical and computational purposes, it is con- 
venient to transform these coordinates to a heliocentric 
ecliptic, mean-of-/r/f/Mch-date system. Thus, two rotation 
matrices must be computed. First, a rotation is made from 
heliocentric equatorial coordinates, mean of 1950.0, to 
heliocentric equatorial, mean of launch date. This is ac- 
complished by means of the matrix A, whose elements are 
given below (as obtained from Ref. 3): 


an = 1 “ 0.000296977^ - O.OOOOOOIST^ 
ai 2 = ~ - 0.022349887 - 0.000006767^ + 0.000002217^ 

a ,3 = - a,, == - 0.009717117 + 0.0000020772 -f 0.000000967^ 
a 22 = 1 - 0.0002497672 - 0.000000157" 

^23 = - 0.0001085972 _ 0.000000037" 

033 = 1 - 0.0000472172 4 . 0.000000027" 


(2) The near-launch-planet phase, in which a launch- 
planet-centered conic is fitted to the heliocentric 
conic for a given launch azimuth and launch site. 

(3) The near-target -planet phase, in which the param- 
eters of the planetocentric conic at the target are 
calculated. 


where 7 is the number of Julian Centuries of 36,525 days 
past the epoch 1950.0. 

The second rotation transforms the coordinates from 
heliocentric equatorial, mean of launch date, to helio- 
centric ecliptic, mean of launch date. This is accomplished 
by means of the matrix E : 


*Coded by W. J. Scholey and R. Y. Roth of the JPL Computer Appli- 
cations and Data Systems Section. 

’Developed by E. Dobies, formerly of JPL, and coded by C. A. Sea- 
feldt, of the JPL Computer Applications and Data Systems Section. 



0 0 
cos € sin e 


" sin e cos € 


(2^2) 
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where € is the mean obliquity of the ecliptic on the launch 
date and sin e and cos € are approximated by the follow- 
ing equations from Ref. 4: 


sin € = sin Co + P{T) cos Co — -- P^{T) sin Co 
cos c = cos Co — P(T) sin Cq — ^ P'^{T) cos Co 


(2^) 


where 


P(T) = - 0.00022711ir - 0.0000000286P 
+ 0.000000008781^ 


(2-4) 


and Co is the mean obliquity of the ecliptic for 1950.0. 

Now, in general, the components of the position vectors 
Rl, Rla, and Rp are found in mean-equinox, ecliptic-of- 
launch-date coordinates by 


Rl 

Rla 

Rp 


EA 


Rl 

RIa 

r; 


(2^) 


The unit vectors Rj^ , RJ may then be found: 

R^' ry 1 — Rr? 


r; 


Rl 


Rl 


Rn 


where 


Rl — |Rl 1 and Rp — |Rjil 


(2-6) 


(2-7) 


The heliocentric central angle i. e., the angle between 
Rl and Rp, is computed from 


Rl • Rp 

^ ^ R R 
ixl rtp 

sin ^ = sgn [(Rl X Rp) • K'] (1 — cos*^ ^)^ 


(2-8) 


where K' is the unit vector normal to the ecliptic plane 
in the direction of the ecliptic north pole. 

The unit vector W normal to the probe's orbit plane, 
is found from the vector equation 


W 


_ Rl X Rp 


(2-9) 


Rl Rp sin 4^ 

Subsequently, the inclination of the orbit to the ecliptic 
plane i may be found from 

cost = W-K' (2-10) 

where 0 t — 7 t/ 2 are the only cases of interest to this 
program. 


Essentially, calculation of the W vector determines the 
planar orientation of the transfer conic. 

2. Application of Lambert's Theorem 
As shown in Section I, the semimajor axis a may be 
found, given the flight time Tp*, Rl, and Rp. One uses 
Lamberts Theorem, which states: 

The transfer time between any two points on an 
ellipse is a function of the sum of the distances 
of each point from the focus, the distance be- 
tween the points, and the semimajor axis of the 
ellipse. 

Or, functionally, 

Tf = T^{Rl + Rp,C,fl) (2-11) 

where C is the chord distance from Rl to Rp, or 

C = 1Rp-Rl| (^12) 

Now, since T y, Rl, Rp, and C are known, Lambert's The- 
orem can be used to solve for a by an iterative process 
described below. 

First, the transfer time Tm of the minimum-heliocentric- 
energy trajc'c tory is computed from Eqs, 2-13 to 2-22. 

The minimum semimajor axis Om for an ellipse can be 
found from 


= (Rl + Rp + C)/4 (2-13) 

and Cmj the eccentricity of the orbit with a minimum semi- 
major axis, is expressed by 


= Ji-^ 

\ Om 


where pm is the semilatus rectum obtained by 
.. _ 2(2fl^ “ Rl) {2am - Rp) 

pm — 


(2-14) 


(2-15) 


The true anomaly at launch VLm is found from Eqs. 2-16 
and 2-17. First, compute the angle <f): 

_ Rl “ ( 1 ~ ^m) 


cos <j) 


Rl 


0 ^ <j) ^ 7T ( 2—16 ) 


where is an angle from the aphelion to Rl. 


It can then be shown from the geometry of the mini- 
mum-energy ellipse (Fig. 2-1) that 


VLm — -JT — </>| 0<^<X 
Olw = 7T-h</)i ?r<^<27r 


(2-17) 
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(a) Centrol angle less than 1 80 deg (b) Central angle greater than 1 80 deg 

Fig. 2-1 . Geometricol configuration of heliocentric conic 


The true anomaly at arrival on the minimum-energy el- 
lipse Vpm can be found from 


^Lm ^ (2—18) 

Now, the mean anomaly at launch in the ellipse can 
be found from 

= Etm — Cm SlTl ELm (^19) 

where Ei^m, the eccentric anomaly, is given by the equations 


cos E^ni - 
and 

sin Eiyri = 


Bm + cos VLm 
1 cos 


yi - sin VLm 
I Cm cos VLm 


0 < Ei,m ^ 27T 


( 2 - 20 ) 


The mean anomaly at arrival Mpm can be obtained in a 
similar manner. Then, the mean-anomaly difference is 


= Mpm — Mi^ (2-21) 

By use of this difference, Tm> the time for a minimum- 
energy trajectory from launch to target, can be found 
from the equation 


Tm = 


AM„ 


27T 


( 2 - 22 ) 


where the units are years and astronomical units (au). 


Next, the semimajor axis a, corresponding to the given 
flight time Tf, is obtained. At this point, there can be 
two types of trajectories, Tf > Tm or Tf < Tm. 


For Tf> Tm, letfli.i = (i + l)a„„ until TF(ai^^) > Tf 


For Tf < Tm, let fli+i = (t -h l)am, until Tjr(ai^i) < Tf 


where t - 1, 2^ ... n, and ai = Om- Note that r^.(ai) is 
the time of flight corresponding to aj, and is calculated 
from Eqs. 2-24 to 2-29 and 2-19 to 2-22; here, however, 
the subscript i is used in place of m. When T(Ci ) < Tf 
< T(ai^i), then Ci < a < Oi^i, and a slope - intercept 
method is used to obtain an approximation on a. Using 
a,^i for a first estimate on a^, and using 


Cj = a — Aa 


Tf ~~ T(flj-i) 

T(a0 - ^ 


Aa - (aj.i - ai) 
subscript ) = 1, 2, ... n 


(2-23) 


for subsequent estimates, the value of a^ is used for a to 
calculate a and ft. 


P = C0S“ 


, r c^ + - Rp^ - 4a (Rp - fl 


2C(2fl - Rl) 

. TRl — Rp cos 4^1 

„ = COS-|_ ^ J 


fO < /? < 7T 
> 

[0 < a < 9T 

(2-24) 


Now, the linear eccentricity X can be found from 
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X2 = 4a2 - 2R^ (2a - Rl) [1 + cos (|a - /J|)] (2~25a) 
X2 ^ 4a^ - 2R^ (2a - Ri) [1 + cos ([a ^ ^])] (2-25b) 

Equation 2-25a is used if Tp > T^, and ^ > tt; or if 
Tp < Tm, and ^ < tt. However, Eq. 2-25b applies if 
Tp > Tm, and ^ < tt; or if Tp < T^, and ^ > tt. 

Then, the true anomaly at launch is found from 
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If Ri, > Rp, Ua and Up are interchanged. The true anomaly 
at arrival Vp is found from 

Op = Vl ^ 0 < Up < 2 7T 


where i is the inclination of the planet's orbital plane to 
the ecliptic and H is the longitude of the planet s ascend- 
ing node. 


and the eccentricity e is found from 


e 


2a 


(2-29) 


Then, T(aj) is recalculated using Eqs. 2-19 to 2-22. 


When \Tf — T{ «i)l < e. a predetermined conver- 
gence criterion, the conic is considered to be the one 
desired, and a{= a»), c,Vl, Vp are the parameters of the 
heliocentric transfer orbit. 


At this point, the transfer orbit is completely deter- 
mined. Other key quantities are calculated as shown in 
the following paragraphs. 

3. Calculation of Velocity Vectors for Probe and Planet 

The heliocentric velocity of the probe V (with the 
subscript L for launch time or p for arrival time ) , can be 
calculated by using the equation 

V = -^[(WXR)cosr + Rsinr] (2-30) 
R 

where 

The path angle F can be found by the expression 


The unit vector P directed toward the perihelion is ex- 
pressed by 


Px = cos o> cos n sin 0 ) sin Q cos i j 

Fy = cos 0) sin Q -h sin <o cos Q cos i > (2-34) 

Pz = sin o> sin i j 

where <o is the argument of perihelion of the planet. 


The unit vector Q right-handed to P and W, is obtained 
from 

Q - W X P 

Then 


cos i; = R' • P 
sin u — • Q 


(2-35) 


where, again, R' is the Sun-planet unit vector. 


Numerical values for the planet elements a, e, i, Q, and 
o) (from Ref. 5) are presented in Table 2-1. 

Finally, the hyperbolic-excess velocity vector at launch 
or arrival V^ can be found from 


sin r 


where 


R 

(1 ~ e^) {2a - R) 


e sin o 


(2-32) 


0<r<|^if 0<u<,r 

-^<r<0 if n <v<2^ 

A 

The velocity vectors Vi and V 2 of the launch and target 
planets, respectively, are found in a similar manner. In 
this calculation, however, the semimajor axis a, the path 
angle F, and the eccentricity of the orbit e all refer to 
the respective planet's orbit around the Sun. The true 
anomaly v of the planet in its orbit at any time can be 
found from Eqs. 2-33 to 2-36. 

The unit vector W normal to the planet’s orbit plane is 
given by 

W = sin n sin i, cos Q sin 1 , cos i (2-33) 


V, = V- Vpiane^ (2-^6) 

4. Calculation of Various Trajectory Parameters 

To assist the trajectory engineer in selecting and de- 
signing interplanetary trajectories, various key trajectory 
parameters are computed. The formulas for these are 
given below. 

The angle y between the hyperbolic-excess velocity 
vector and the planet’s orbital plane can be found at 
launch and arrival by 

y ^ y ^ 2 ^ 2—37 ) 

The right ascension 0 and the declination 4> of the 
asymptote (launch or arrival) can be found by use of 
the expressions 


sin ^ = Sz “ < 4> < ^ (2-38) 


If 
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Table 2-1. Mean planet elements 


Planet 

Semimajor oxis a, 
au 

Eccentricity e 

Inclination to ecliptic ir 
deg 

longitude of ascending node Q, 
deg 

Argument of perihelion u), 
deg 

Mercury 

Venus 

Earth 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

Pluto 

0.387098 
0.723331 
1 .000000 
1.523679 
5.2027 
9.546 
19.20 
30.09 
39.5 

0.205625 + 0.0000207“ 
0.006793-0.0000507 
0.016729-0.0000427 
0.093357 + 0.0000947 
0.048417 + 0.0001647 
0.055720-0.0003457 
0.0471 +0.00027 

0.00872 +0.000047 
0.247 

7.003819 + 0.001757 
3.394264 + 0.001257 
0 

1.849986-0.0006397 
1.305875-0.0056947 
2.490583-0.0038897 
0.772792 + 0.0006397 
1.774486-0.009537 
17.140000-0.005567 

47.737778 + 1.185007 
76.236389 + 0.905567 
0 

49.173611 +0.773897 
99.94861 1 +1.010567 
1 13.226806 + 0.873067 
73.726528+0.498617 
131.230833 + 1.098897 
109.633750 + 1.358067 

28.937778+0.369447 
54,619305 + 0,501397 
102.078056 + 1.716677 
285.965000 + 1.066677 
273.577222+0.599447 
338.850694 + 1.085287 
96.129028 + 1.1 12507 
272.935934-0.432227 
1 13.860694+0.030837 

^Measured ii 

n Julian Centuriei from 19 

50.0 


and 


cos 0 


sin 0 


__h_\ 

ys% + s?, / 

Sy ( 

\/S^ + Sy I 


0 < 0 < 27t 


(2-39) 


The angle between the launch hyperbolic-excess 
velocity and the Sun-launch-planet line at launch is de- 
fined by 


cos u ^ 0<U<- (2-43) 

y hL 


where 

S (S„Sy, S^) (2-40) 

V h 

Here, £ is the rotation matrix given by Eq. 2-2; S is a 
unit vector in the direction of the outgoing asymptote 
when \h is calculated at launch, and in the direction of 
the incoming asymptote when \h is calculated at the 
target. 


There are six other angles • Cpj the Sun-probe-target 
angle; the Sun-target-Earth angle; fr, the probe- 
tar get-Canop us angle; rjsy the supplement of the angle 
betw^een the projection of the target-Sun vector on the 
R-T plane and the T direction; the supplement of the 
angle between the projection of the target-Earth vector 
on the R-T plane and the T direction; the supplement 
of the angle between the projection of the target-Canopus 
vector on the R-T plane and the T direction. Here, T is 
a unit vector lying in the ecliptic plane, given by 


Other quantities which are of interest for the helio- 
centric phase are given in Eqs. 2-41 to 2-48, below. 

The communication distance at arrival Rr is expressed 
by 


Rc — |Rc| 

Rc “ Rp Rl4 


and 

(2-41) 


The arrival ( or departure ) angle a is obtained from 


cos a 


^planet * 
^planet 


0 ^ (X ^ 7T 


(^2) 


S, XK' 
[SpXK'l 


R Sp X T 


(2-44) 


cos fp = - Sp • RJ 1 

cos = — Sp • RJ, > 0 < f < 7T 

cos ic ~ Sp*C j 

(2-45) 
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where C is a unit vector to the star Canopus, and is 
obtained by normalizing Rc from Eq. 2-41. 


sin 1 ), — 
cos rjs — 
sin r}B — 

COS rjE — 

sin TfQ — 


RR^ 

sin^:p 

TR‘ 

sinCp 

RR^ 

sin^:^ 

T*Ri 

sin (n 

RC 

sinfo 



COS Tjc — 


T»C 

sin fc / 


(2-46) 


The angle ap between the projection of the incoming 
asymptote on the target planet s orbital plane and the 
target-planet-Sun line at arrival time is defined by 

cos <Tp — —Rj, • Spr ) 

> ““R' ^ cTp < IT (2-47) 

sincTp == - Spr* (W2 X R^) j 


where Spr, the projection of Sp on the target planet's 
orbital plane, can be found by 


. _ Sp- W2(Sp-W2) 
|Sp- W2(Sp-W2)l 


(2-48) 


Here, is a unit vector normal to the target planet's 
orbital plane. 


6. The Planefocentric-Conic Trajectories 

The trajectory near the launch planet is found by 
application of several conditions: 


(4 ) The parking-orbit altitude plus the launch planet's 
radius is equal to the perifocal distance of the 
escape hyperbola. This is a good practical approxi- 
mation, since it is most eflBcient to inject the space- 
craft into the escape hyperbola near the perifocus. 

(5) The launch planet is assumed to be spherical in 
shape. 

Given these conditions, the following formulae are used 
to compute the parameters of the near-launch-planet 
trajectory. 


The eccentricity e of the launch-planet conic can be 
computed from 


e = 1 + 


RpCs 

GMl 


(2-49) 


where Rp is the perifocal distance of the near-launch- 
planet conic, and C 3 is the vis viva energy, defined as 


yz 

^ hL 


(2-50) 


and GMl is the universal gravitational constant times the 
mass of the launch planet. The radius to injection R is 
found from 


R - 


P 

1 + e cos V 


(2-51) 


where v is the true anomaly of injection on the launch- 
planet escape hyperbola, and p is the semilatus rectum 
given by 


-GM^(1 -e^) 
C 3 


(2-52) 


The path angle at injection F is found from 


cos r = 


\/pGMl 

VR 


o<r< ^ ,ifo<t)<7r 

-^<r<0,if7r<u<2T 

(2-53) 


(1) The injection energy C 3 and the direction S of 
the outgoing asymptote of the escape hyperbola 
are obtained from the solution of the heliocentric- 
transfer problem as given in Section II-A above. 

( 2 ) The vehicle is launched from a given launch site, 
specified by its latitude and longitude, into a 
low-altitude circular parking orbit. 

(3) After coasting in the parking orbit for a time tc, 
as determined below, the final stage ignites and 
propels the spacecraft to the final injection energy. 


where V, the injection speed, is given by 

V = yjcs + (2-54) 

The vector W, perpendicular to the plane of the conic, 
can be found by the solution of the two vector equations 

W • S = 0 ) + WySr + WzS^ = 0 

W • W = 1 f W^;i + + W% = I 

(2-55) 
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where 


_ -(WySy + WzSz) 

\y ^ ^ 


(2-56) 


and 


Wy 


W^zSySz 


Sj: 

S| + S^ 


VI - si - wi 

(2-57) 


It should be recalled that S is a unit vector in the direc- 
tion of the outgoing asymptote of the escape hyperbola. 
Thus, a condition is set on Wz: 


W| < 1 - Si (^8) 


From the geometry (Fig. 2-2) of the launch azimuth 
and the launch latitude ^>l, it can be shown that 

Wz = cos sin ( 2-59 ) 

From Eqs. 2-^8 and 2-59, the W restrictions can be 
calculated as a restriction on Xl, since and S are fixed: 

sin^S^<Vp- ( 2 - 60 ) 

COS^^L 

A vector B, orthogonal to S and W, is calculated: 

B = SXW (2-61) 

Now, let P be a unit vector in the direction of the peri- 
focus as shown in Fig. 2-3: 

P S cos t)^ -h B sin Vg ( 2-62 ) 



Fig. 2-2. Trajectory-plane-launch-site geometry 



Fig. 2-3. Trajectory-plane geometry 


Then Q is the unit vector at right angles to P; 

Q = S sin u, ~ B cos Vg ( 2-63 ) 

The true anomaly of the outgoing asymptote S is 
given by 

cos u, = — ^ 0 < Vg < w (2-64) 


The right asc'f^nsion 0^ of the launcher can be found from 


cos 0L = 


sin 0L = 


Wx sin sin + Wy cos Xl 

W| ~ 1 

Wy sin sin Xl ^ Wx cos Xl 
W| - 1 


0 < 0Z. < 27T 

(2-65) 


where is the latitude of the launcher. A unit vector R ) 
in the direction of the launcher can be calculated in the 
vernal-equinox equatorial system by 


Rl — cos cos 0L, cos sin 0 l, sin 


(2-66) 

The angle ^ between the launcher and the perifocus of 
the conic is given by 

cos = RV P ) 


sin 




0 < ^> < 27T 


(2-67) 


The angle between launch and injection is defined by 


^/ = 27r-$ + t; (2-68) 

where, again, v is the true anomaly at injection into the 
escape hyperbola. 
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A unit vector toward injection is seen in Fig. 2-4 to be where 


R^ = P cos u -f Q sin t; = ( Rx, Ry, 

Rz) 

(2-69) 

it 

is the time of first burn ( into the parking orbit ) 

The injection latitude is given by 



^2 

is the time of second burn ( into the escape 
hyperbola ) 

sin Rz — ^ < 

7T 

(2-70) 


is the angle of first burn 


2 

<J>2 

is the angle of final burn 

The right ascension of injection 0 can 

be found from 


is the inverse parking-orbit rate = 14.689 sec/deg 


for a 100-nm Earth parking orbit 




Fig. 2-5. Powered trajectory modified by coasting 


The azimuth at injection X can be computed from 
Sz — ks sin ^ 


cos X = 


ki cos 4> 


0 21 <C 7T 


where 


(2-75) 


, Ve^TTX . 1 

ka = Sin u cos V 

e e 


j - 1 , 1 . 

COS V -t — sin V 

e e 


(2-76) 


The unit vector to final-stage ignition can be found in 
a manner similar to that used for R: 


Fig. 2-4. Ascent-trajectory geometry 


Rz = P cos ( u — ^ 2 ) + Q sin ( u — ^ 2 ) (2-77 ) 
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The latitude <E>2 of final-stage ignition is computed from 
sm^2 = Rz2 — (2-78) 

The right ascension 02 of final-stage ignition is given by 


cos 02 


sin 02 


Rx2 \ 

VRL + Rk / 

R ) 0 < 02 < 27T 




(2-79) 


The injection time T is calculated from 

r = + (2-83) 


C. Differential Corrections 

As outlined in Section I-C, the calculation of difiFeren- 
tial corrections, or partial derivatives relating variations 
in the impact parameter B and flight time Tf to variations 
in the hypei bolic-excess velocity vector V/,/, at launch, are 
performed by a numerical differencing technique.® 


and the longitude ^2 of final stage ignition is defined by 

^2 — ^ "b 02 — 0 <ot2 (2—80) 

The angle Av between injection and the outgoing asymp- 
tote is expressed as 

Av = Vs — V 

The launch time Tz, can be found by 

O<(0i - - GHA)<2»- 

(2-81) 

where GHA is the Greenwich Hour Angle at 0^ UT of the 
launch day as shown in Fig, 2-6, and is obtained from 

GHA - 100? 07554260 

+ 0?9856473460 T, 0 < GHA < 2n (2-82) 
+ 2?9015 X 10-'^ Tl 

Here, Td = days past 0* January 1, 1950. 



The basic idea in this technique is to compute a varied 
or perturbed trajectory and then difference it with the 
reference case. A small variation A\ hi in the hyperbolic- 
excess velocity vector is equivalent to a small variation 
A\ L in the launch heliocentric-velocity vector. Letting 
primed quantities denote variables on the perturbed 
trajectory, the launch heliocentric velocity on this trajec- 
tory is, then, 

V' + AV,, (2-84) 

where 


aV„ = (C:)^ A<^s[ — sin ^scos0s, — sin^.sjsin0.s, cos<^.sd, 
( C i A0X [ — cos sin 0s, cos cos 0,?, 0] , 

- [ cos cos 0s, cos sin 0s, sin ^s] 

2 (C;d^ 


where A0s are small angular variations (0.2 deg), 
and the energy variation is AC3 = 0.005 C3. 


The semi major axis a' is obtained from 

R. 


a' =■ 


GMs 

The radial rate R,' is given by 

• V'*Rr 

R' = t 

^ Rl 


(2-85) 


( 2 - 86 ) 


The semilatus rectum p' and eccentricity are computed 
from 




Rl (vr - R^ 

GMs 


(2-87) 


€f = 



( 2 - 88 ) 


Pig. 2-6. Relation between longitude^ right ascension, and 
time in equatorial plane 


This method was developed by William Kizner, Research Specialist, 
Systems Analysis Section, Jet Propulsion Laboratory. 
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The eccentric anomaly at launch E' is expressed by 


sin E' = 

u 


Rl r1 


e' (a'GMs)^ 


The mean anomaly at launch is obtained from 
= E' — e' sin E' 

The mean orbital rate n' is given by 


(2-89) 


(2-90) 


, „ (GM,)^ 


n" ~ 


rj^3/2 


The mean anomaly at the target M' is stated as 
M' = n' Tj. -f M' 


(2-91) 


(2-92) 


The eccentric anomaly at the target E'p is obtained from 
the expansion 


e;=e,+ 


( I ^ 

yi — e' cos Ep/ 


am 


. cos E' —e' 
cos V' = p 

1 — e' cos Ep 

. , (1-O'^^sinE' 

sin v' = 2 1 L 

1 — e' cos E' 


(2-96) 


The heliocentric central angle is given by 


p 


(2-97) 


(2-98) 


The angular momentum h' is expressed as 
h' = Rt X V' 

The heliocentric position vector at the target is given by 
RJ = R; (^ cos <!r' + sin ❖') ( 2-99 ) 

where 

R; =fl'(l - e'cosEp (2-100) 

A vector in the direction of perihelion with magnitude o' 
is computed from 


if 


, 1 r (Sg'sinEp)^ — (1 

eL (1- 


If e' sin Ep 

AM= 

~2 L(1 - e'cosEpY 

— e' cos Ep) (e' cos Ep) ' 

AM^ 

e' cos Ep)® 


(2-93) 


cos Ep ^ 0 


However, 

e; = e. 


^ e cos Ep — 1 + -\/Je cos Ep — 1)^ + (2e sin Ep) AM 

e sin Ep 


if 


(2-94) 


cos Ep < 0 

where 

AM = Mp — (Ep — e' sin Ep) 

The true anomalies at launch and at the target, u' and c' 
are found from 


cos o' 


ry-RL 

e'Ri 


(2-95) 


0 < u' < 7T if R ' is positive 

* 

^ < o' < 2^ if Rl is negative 


V'i X h' 
GMs 


- If (^101) 


The heliocentric velocity at the target is defined by 


v;- 


= 


(t-') 


(2-102) 


The hyperbolic-excess velocity at the target is expressed 
by 

Vip = v; -V, (2-103) 

The diflFerence between the heliocentric position vectors 
on the perturbed and reference trajectories is given by 

aR;=R'-Rp (2-104) 

The impact parameter R is computed from 

„ (aR' •¥')¥' , 

B = — f hp/ hp ^ ^IJ/ 

VZ 

The flight-time error is stated as 

aR: • v;„ 

aTe = -- V'» (^105) 

^ hp 


The partial derivatives are formed by dividing A0s, 
A4>s, and AC3 into the miss components B • T, B • R, and 
flight-time error ATp. In addition to the component 
partials, the quantity dB/dQi is defined by 


dB 

dQi 


"/aB-TV /aB*Ry“|^ 

A ) \ . 


(2-106) 
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The three partials, dB/d&a, 0B/04>s, 0B/0C , are im- 
portant measures of the error sensitivity of a trajectory. 


where an is tlie astronomical-unit-to-kilometer conversion 
factor 


The eflFect of uncertainty in the knowledge of the 
astronomical-unit— to— kilometer conversion factor on target 
miss and flight time may be determined by the following 
formulae: 


whence 


and 


aB*T -2C3 dB T 

dau au 9C3 

aB'R „ -2C, 6iB-R 

dau an dCz 


SB _ 2C, dB 
dau au dCs 


dTp _ -2C3 dip 
dau ^ au 6C3 


(2r-107) 


(2-108) 


(2-109) 


The effect of solar-radiation pressure acting on the 
probe may also be evaluated as follows: In Eq. 2-84, let 
AV/i/, “ 0, but in Eqs. 2—85, 2-87, 2-89, 2r-91, and 2-101, 
vary GMs by adding an increment AGMs. This procedure 
gives rise to a varied trajectory from which the impact 
parameter B and flight-time error aTp may be obtained. 
The partials oB/dGMi^ and cTp/cGMs may then be cal- 
culated. Since the acceleration caused by solar-radiation 
pressure acts in a direction opposite to the gravitational 
attraction of the Sun, radiation pressure has the effect of 
decreasing the Sun's gravitational attraction or 
decreasing GAis. A decrease, AGM.s- = —2.4 X 10^ km^/ 
sec-, corresponds to the solar-radiation pressure acting on 
a 300-kg spacecraft having a perfectly reflecting area of 
3.6 m^ Thus the miss, always being a positive number, is 
obtained by AB,p 2.4 X 10" dB/dGMs, and the cor- 
responding f light-time error is = — 2.4 X 10*" 

dTi/dGM;<, which is sign-sensitive. 
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III. DISCUSSION AND EXPLANATION OF RESULTS 


A. Introduction 

1. Trajectory Computations 

Minimum (geocentric) energy trajectories from Earth 
to Mercury for the period 1964-1975 and from Earth to 
Jupiter for the period 1968-1973 were computed on the 
7094 digital computer by numerically evaluating the 
analytical model explained in Sections I and II. This 
computation resulted in curves of C , (vis viva geocentric 
energy) vs launch date for these trajectories. 

For each target planet, several 75-day launch intervals 
were analyzed in greater detail: curves of 19-21 key tra- 
jectory parameters vs launch date for various ri.s viva 
geocentric energies (C.) were automatically plotted. 

2. Trajectory Analysis 

Careful analysis of the results shows that as many as 
four ballistic flight paths to the target planet exist per 
launch date for a given r/.v viva geocentric energy (C;), 
assuming trips of less than one revolution around the 
Sun. A more extensive analysis might reveal more than 
four flight paths. Actually, feasible launchings can occur 
only for small time intervals (1-3 months), when the rela- 
tive i^ositions of the Earth and target planet are such 
that the velocity requirements for ballistic transfers can 
be reasonably achieved by modern boost vehicles. These 
intervals occur once during each synodic period of the 
planet. A synodic period is the time interval required for 
tlie Earth and target plaiu't to attain successive identical 
relative angular relationships in heliocentric longitude, 
e.g., the time between Earth and target planet align- 
ments (116 days for Mercury; 390 days for Jupiter). 

Figures 3—1 to 3—6 show the minimum injection energy 
recpiired for transit to Mercury vs launch date. The dates 
of Earth-Mercury alignment are indicated by arrows 
(t). Favorable launch opportunities (i.e., relative minima 
in injection energy) and planet alignments occur about 
every 116 days. Two families of conic trajectories exist 
(Types I and II), and will be describe^d in detail in Sec- 
tion B. After three synodic periods (348 days), Mercury 
and Earth again assume approximately the same space- 
fixed position. This condition is also reflected in Figs. 
.3-1 to .3-3, in which the minimum energy curve for 
Type I trajectories has very pronounced minima at .348- 
day intervals. Type II trajectories exhibit this periodic 
behavior also, although the curves are more complicated. 


Within the time period considered, the minimum injec- 
tion energy reaches an absolute minimum on November 
23, 1967, for Type I trajectories, and on November 3, 
1968, for Type II trajectories. Only the launch oppor- 
tunities centered around these dates are discussed in 
detail in Sections IV and V. The results are fairly repre- 
sentative for all other launch opportunities. 

If the synodic period of the Earth and target planet 
were a rational fraction of a year, and if their orbits 
remained unchanged, the minimum energy cairves would 
be exactly periodic functions of time. 

Figures .3-7 and .3-8 show the minimum injeetion energy 
for transits to Jupiter vs launch date for Type I and 
Type II trajectories. The dates of Earth-Jupiter alignment 
are indicated by arrows (|). 

Favorable launch opportunities, planet oppositions, 
and approximately the same space-fixed geometry of 
Earth and Jupiter occur every synodic period (390 days). 
For launch opportunities, the absolute minimum injec- 
tion energy occurs on December .30, 1969, for Type II 
transfers, and on January 3, 1970, for Type I trajectories. 
The launch opportunities centered around these dates 
are discussed in detail in Section VII. In addition, be- 
cause of some interesting features that were discovered, 
the launch opportunity extending from November 1968 
to January 1969 is treated in detail in Section VI. 

B. Classification of Trajectories 

1. Type I and Type II 

In observing the variation of any trajectory parameter 
vs launch date for fixed geocentric energies C.^, it is noted 
that two separate groups of closed contours, rather than 
one, are described (see Fig. .3-9). These two sets of energy 
contours are designated as Type I and Type II trajec- 
tories, where Type I trajectorit's are defined as having 
heliocentric-transfer angles less than 180 deg, and Type II 
trajectories are those having heliocentric-transfer angles 
greater than 180 deg. A Type I trajectory traverses less 
than half-way around the Sun from launch to planet 
encounter, whereas a Type II trajectory would traverse 
more than half-way, but less than one full revolution, 
around the Sun. For a given launch day and energy, then. 
Type II trajectories require greater flight times than do 
Type I. The existence of these two sets of energy contours 
can be attributed to the fact that the orbits of Earth and 
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. Mercury 1964-1968: Minimum injection energy vs launch date, Type I 
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Fig. 3—2. Mercury 1968—1972: Minimum injection energy vs launch date, Type I 
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Fig. 3—4. Mercury 1964-1968: Minimum injection energy vs launch date, Type II 
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Fig. 3-5. Mercury 1968-1972: Minimum injection energy vs launch date, Type II 
























LAUNCH DATE 




JPL TECHNICAL REPORT NO. 32-77 


the destination planet are not coplanar. If the Earth and 
tar^(‘t planet had coplanar circular orbits, only a single 
group of (‘iKTgy contours would exist. If the orbits of 
Earth and the target planet were circular as well, this 
singl(‘ group of energy contours would include the well- 
known Ilohmann’s tratisfer orbit. The Hohmann transfer 
orbit is tangential to the orbits of Earth and target planet, 
lias a lu'liocentric transfer angle of 180 deg, and requires 
the least vis viva geoc(‘ntric energy (C ;) among all transfer 
trajectories. TIk' Hohmann minimum-energy trajectory 
would iH'main constant for eac h synodic period, retaining 
tlu‘ same flight time and vis viva geocentric energy. The 
fact that tile orbits of Mercury and Jupiter are slightly 
incliiu‘d with respc'ct to the I'cliptic (7 and 1.3 deg, re- 
spectively) strongly influences the* transfer trajectories 
and causes the single group of trajectory contours to split 
into two families, tlius necevssitating the dc^signation of 
Typ(^ I and Type* 11 trajectories. 

In Fig. 3-9 note the small family of closed contours 
ccaitered around January 1, 1969. This group has been 
designated Type I- A. Type' I-A trajectories exist for 
Venus also. They are' neit attractive because their flight 
times are longer for a given vis viva geocentric energy. 

2. Class I and Class II 

As Fig. 3-9 sheiws, twei trajectories usually exist for a 
given type, energy, and launch date. Trajectories correv 
sponding to the lowec portions of the energy contours 
are (h'signated as Cflass I and are' connected by solid 
lines and corrc'spond to shorter flight times and smaller 
heliocentric central angle. Trajtctories corresponding to 
tlu‘ upp('r part of the eneegy contours are designated as 
Class II and are* eomucteel by dashed lines and corre- 
spond to longeT flight timers anel larger heliocentric angle. 
For a given launch date and ge'ocentric energy, the Type 
I-Class I trajectory has th(' shortest flight time and the 
Typ(‘ 11-Class II trajectory has the longest flight time. 

For example, Fig. 3-9 shows that for the launch date 
Noveml)er 26, 1969, and C , - 120 kmVsec^ the Type I- 
Class I trajectory has a flight time of .520 days, whereas 
the Type ll-Class II trajectory has a flight time of 1460 
days. The flight time for the Type 1-Class II trajectory is 
880 days, and for the Type 1 1-Class I trajectory it is 982 
days. 

3. Minimum Energy Trajectories 

a, Minitnum-energij loci. Further study of Fig. .3-9 
shows that the geocentric energy C t has a relative mini- 
mum of C > 85.617 kmVsec“ on December 4, 1968, 

where the flight time is 8.50 days. By moving away from 


this point in any direction, increases. Type I trajec- 
tories corresponding to - 87 km-/sec- only exist from 
November 29, 1968 (flight time 780 days) to December 
10, 1968 (flight time 950 days). The geocentric energy 
required is smaller inside the contour and larger outside 
it. Consequently, the minimum possible energy on No- 
vember 29 and December 10, 1968, for Type I trajectories 
is 87 krn‘/sec“. A curve of minimum geocentric energy 
vs launch date may be constructed for Type I and Type II 
trajectories by constructing tlie vertical tangents to the 
contours of a given parameter such as launch time. By 
connecting these points, a locus of minimum energy 
(Fig. .3-9) is obtained which separates Class I and Class II 
trajectories within a given type. 

Figure .3-10 was constructed in this way from Fig. 3-9 
and shows minimum geocentric energy vs launch date 
from November 18, 1968, to January 7, 1969. The curve 
corresponding to Type II trajectories was terminated on 
December 28, 1968, because flight times corresponding 
to later launch dates are in excess of four years. Two 
relative minima, C ^ ~ 85.617 kmVsec- on December 4, 
1968, and 89.410 kmVsec ' on January 1, 1969, exist 
for Type I trajectories because of the noncoplanar and 
noncircular nature of the orbits of the Earth and Jupiter. 
The relative minimum for Type II trajectories is = 
77.8.32 kmVsec- and occurs on December 13, 1968. 

b. Absolute minimum-energy trajectories. The abso- 
lute minimum-energy trajectory for each type is defined 
as the trajectory with the least geocentric energy within 
each synodic period of the planet. Such trajectories are 
useful as first approximations of launch dates, communi- 
cation distances at planet encounter, and flight times for 
missions within a given launch opportunity. Table ,3-1 
shows the characteristics of absolute minimum -energy 
trajectories to Mercury and Jupiter. These trajectories 
occur once every synodic period of the target planet, as 
explained in Section III. A. 2. For Mercury, a sharp de- 
crease in the absolute minimum energy occurs every third 
synodic period (348 days), as can be seen from Figs. .3-1 
to .3-7. Table .3-1 shows that the lowest possible geo- 
centric energy to Mercury is 41.2 kmVsec- (launch date 
October 31, 1968, Type II); to Jupiter it is 75.3 kmVsec- 
(launch date January 3, 1970, Type I). 

Very recently Cutting and Sturms (Ref. 6) showed that 
it is possible to reduce the energy requirements for a 
Mercury mission by 70% by an encounter with Venus 
before arrival at Mercury. The following is a typical 
trajectory: Launch date August 14, 1970, injection energy 
13 kmVsec-, arrival date at Venus, November 26, 1970, 


29 






JPL TECHNICAL REPORT NO. 32-77 


Table 3-1. Characteristics of minimum-energy-transfer 


Planet 

Trajectory 

type 

Launch 

date 

Flight 

time, 

days 

Geocentric 

injection 

energy,' 

m‘Vs‘^ X 10^ 

Heliocentric 

central 

angle,' 

deg 

Sun>planet 
distance at 
arrival, 
km 

Earth-planet 
distance at 
arrival, 
10^ km 

Celestial 
latitude 
of planet 
at arrival, 
deg 

Mercury 

1 

2/1 1/64 

95 

0.593 

128.4 

69.3 

102 

-4.69 


1 

6/14/64 

95 

0.895 

153.6 

46.4 

135 

1.04 


1 

10/6/64 

86 

1.006 

142.6 

53.2 

126 

6.70 


II 

3/1/64 

105 

0.926 

215.7 

50.5 

178 

-3.79 


!l 

6/18/64 

105 

1.053 

234.2 

50.9 

188 

6.99 


II 

9/24/64 

116 

0.497 

221.5 

66.8 

176 

0.70 


1 

1/21/65 

97 

0.517 

1 35.6 

69.8 

106 

- 3.25 


1 

5/29/65 

97 

0.847 

163.7 

46.8 

142 

0.43 


1 

9/19/65 

86 

1.074 

139.9 

50.2 

125 

7.04 


1 

12/31/65 

101 

0.467 

146.6 

69.5 

1 12 

-2.09 


II 

2/8/65 

109 

0.935 

218.2 

53.6 

178 

-5.52 


II 

5/27/65 

107 

1.085 

216.1 

46.6 

174 

5.60 


II 

9/10/65 

1 15 

0.561 

229.8 

65.7 

184 

1.40 


II 

1 1 /16/65 

150 

0.861 

203.6 

69.8 

121 

-3.32 


1 

5/15/66 

97 

0.829 

172.2 

47.1 

150 

-0.24 


1 

9/3/66 

66 

1.095 

139.6 

48.2 

126 

6.65 


! 

12/11/66 

104 

0.434 

1 56.4 

68.5 

119 

-0.83 


II 

5/6/66 

107 

0.893 

187.0 

46.7 

154 

0.50 


II 

8/25/66 

115 

0.636 

2 36.4 

64.0 

189 

2.40 


II 

1 1/8/66 

1 35 

0.646 

184.1 

67.7 

1 15 

-0.1 1 


1 

4/26/67 

92 

0.897 

133.7 

55.1 

1 20 

-6.07 


1 

8/17/67 

87 

1.077 

1 39.2 

46.8 

126 

5.66 


1 

1 1 /23/67 

107 

0.41 2 

169.6 

67.8 

1 30 

-0.16 


II 

Al25l67 

104 

0.842 

192.6 

47.1 

161 

-0.20 


II 

%l\2l67 

1 1 2 

0.722 

239.4 

62.0 

192 

3.55 


II 

y\l7l67 

1 23 

0.470 

185.7 

67.8 

130 

-0.16 


1 

4/5/68 

91 

0.834 

122.3 

61.6 

108 

-6.98 


1 

7/28/68 

90 

1.039 

139.9 

46.1 

126 

4.12 


1 

1 1/12/68 

103 

0.449 

177.4 

67.5 

143 

0.10 


II 

4/12/68 

102 

0.858 

199.3 

47.5 

168 

-0.85 


II 

7f7%j6^ 

1 10 

0.822 

242.8 

59.6 

195 

4.61 


III 

1 1 /I /68 

1 15 

0.410 

191.3 

67.9 

146 

-0.24 


1 

3/15/69 

92 

0.738 

1 20.7 

66.2 

102 

-6.42 


1 

7/1 1/69 

92 

0.987 

144.6 

46.0 

1 29 

2.91 


1 

11/1/69 

90 

0.778 

153.2 

60.6 

131 

4.15 


II 

3/27/69 

103 

0.891 

204.5 

48.6 

171 

-2.15 


II 

7l\2l69 

108 

0.921 

241.9 

56.4 

193 

5.86 


II 

10/18/69 

1 15 

0.428 

203.0 

67.5 

159 

0.10 


1 

7l22f70 

94 

0.644 

1 24.7 

68.6 

101 

-5.37 


1 

6/24/70 

94 

0.929 

1 49.4 

46.3 

132 

1.64 


1 

10/16/70 

87 

0.940 

1 45.9 

55.7 

127 

6.05 


II 

3/10/70 

105 

0.910 

21 1 .4 

49.9 

175 

-3.24 


II 

6/28/70 

106 

1.014 

239.2 

53.1 

191 

6.75 


II 

10/3/70 

1 16 

0.468 

215.6 

67.1 

171 

0.38 


1 

2/1/71 

96 

0.562 

1 29.9 

69.7 

102 

-3.93 


1 

6/7/71 

96 

0.873 

154.4 

46.9 

135 

0.17 


1 

9/29/7! 

86 

1.041 

139.9 

51.5 

124 

7.00 


II 

2/23/71 

105 

0.927 

217.0 

51.4 

179 

-4.25 


II 

6/9/71 

106 

1.075 

226.6 

48.6 

182 

6.73 


II 

9/19/71 

I 1 5 

0.524 

223.9 

66.1 

179 

1.10 


1 

1/13/72 

98 

0.497 

138.8 

69.7 

108 

-2.77 


1 

5/22/77 

97 

0.836 

163.7 

47.3 

143 

-0.46 


1 

9/11/72 

86 

1.090 

136.0 

48.7 

123 

6.84 


1 

12/22/72 

102 

0.451 

149.8 

69.0 

114 

-1.43 


II 

1/27/72 

1 14 

0.933 

221.8 



55.5 

- 

177 

-6.16 
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Table 3-1 (Cont^d) 


Planet 

Trajectory 

type 

Launch 

date 

Flight 

time, 

days 

Geocentric 
injection 
energy/' 
mVs2 X 

Heliocentric 

central 

angle, 

deg 

Sun-planet 
distance at 
arrival, 
10^ km 

Earth-planet 
distance at 
arrivol, 
10^ km 

Celestial 
latitude 
of planet 
□t arrival, 
deg 

Mercury 

II 

5/7/72 

113 

0.961 

184.2 

46.8 

147 

0.30 


II 

9/4/72 

114 

0.592 

232.0 

65.1 

187 

1.83 


II 

11/1 2/72 

146 

0.782 

201.7 

69.7 

123 

-2.79 


1 

5/13/73 

93 

0.858 

173.3 

47.2 

154 

- .29 


1 

8/26/73 

87 

1.098 

140.7 

47.6 

127 

6.40 


1 

12/3/73 

105 

0.423 

161.3 

68.1 

1 23 

-0.45 


II 

5/4/73 

103 

0.858 

188.1 

46.8 

1 58 

0.46 


II 

S/20/73 

1 13 

0.672 

237.4 

63.2 

191 

2.85 


II 

11/8/73 

129 

0.555 

183.6 

67.7 

121 

-0.10 


1 

4/18/74 

91 

0.874 

128.2 

57.8 

115 

-6.66 


1 

S/9/74 

88 

1.067 

139.9 

46.4 

127 

5.1 8 


1 

1 1 /I 6/74 

108 

0.412 

176.3 

67.7 

137 

-0.09 


II 

4/22/74 

101 

0.842 

194.6 

47.2 

165 

-0.20 


II 

8/6/74 

1 1 1 

0.762 

241.2 

61.2 

194 

3.92 


II 

1 1/8/74 

1 16 

0.419 

184.3 

67.7 

137 

- 0,09 


1 

3/27/75 

92 

0.798 

1 21.4 

63.9 

105 

-6.92 


1 

7/22/75 

90 

1.023 

1 39.4 

46.1 

126 

3.54 


1 

n /I 3/75 

97 

0.555 

177.0 

67.4 

150 

0.18 


II 

4/5/75 

103 

0,873 

200.6 

48.1 

168 

- 1.56 


II 

7/24/75 

108 

0.862 

243.2 

58.8 

196 

4.92 


1! 

10/28/75 

114 

0,414 

195.9 

67.8 

152 

- 0.17 

Jupiter 

1 

1/2/70 

985 

0.752 

178.8 

779.0 

742 

0.00 


II 

12/31/69 

994 

0.753 

1 81.4 

778.7 

758 

-0.01 


1 

12/4/68 

852 

0.856 

165.8 

804.4 

703 

0.88 


lA 

12/31/68 

1349 

0.884 

1 80.4 

779.2 

735 

0.01 


II 

12/13/68 

1 277 

0.778 

1 91 .3 

784.1 

636 

0.1 8 


1 

1/31/71 

808 

0.777 

1 67.8 

767.4 

783 

-0.41 


II 

2/6/71 

1180 

0.833 

194.8 

750.0 

816 

- 1.02 


1 

3/6/72 

744 

0.812 

162.0 

751.6 

882 

- 0.96 


II 

3/31/72 

1396 

0.857 

198.8 

742.1 

771 

- 1.23 


1 

4/12/73 

721 

0.841 

159-9 

741.7 

889 

-1.29 


II 

5/1 1/73 

1404 

0,839 

196.3 

754.4 

812 

- 0.73 

=* Twite the t< 
^The angle » 
^ Meotured ff 

jtol energy per unit mosi, or yis viVa integral. 

ubtended at the Sun between the Sun-Earth line at launch c 

om the ecliptic. 

nd the Sun-planet line at arrival time. 





Venus perigee altitude 2676 km. Mercury arrival date 
January 30, 1971, hyperbolic excess velocity at Mercury 
12.3 km/sec. 

For transfers to Mercury, the geocentric energy is 
smallest if (1) Mercury is at aphelion (66.8 X 10’’ km) at 
encounter, (2) Mercury is in the ecliptic at encounter, 
(3) the heliocentric central angle is 180 deg, and (4) Earth 
is at perihelion at launch. Owing to the large inclination 
(i == 7 deg) and eccentricity (c -- 0.2) of Mercury's orbit, 
large variations in absolute minimum geocentric energy 
occur, as Table 3-1 shows: The trajectory (launch date, 
October 31, 1968) corresponding to the lowest geocentric 
energy (C.^ = 41.2 kmVsec") has a heliocentric central 


angle of HCA 192.3 deg, the Sun-Planet distance at 
arrival is S[M^ - 67.9 X 10'' km, and tlu‘ celestial latitmh^ 
at arrival is PAP -- 0.46 deg. The trajectory with launch 
date Kovi^mhcv 24, 1967, has almost th(‘ same encagy 
(41.3 km-/see“), even though the lu4iocentrie transt(T 
anglt' is smallcT by 24 deg (HCA 168.6 dc'g). 1 h(^ 
trajectory launehc“d on May 15, 1966, shows the (*ffect of 
arrival near Mercury perihelion (SPP -- 47.6 X 10’ km). 
The geocentric energy is nearly doubled (C ; 84.8 

km-/sec“). It is difficult to find a trajectory in Table 3-1 
which would show only the effect of reaching Mercury 
out of the ecliptic. However, several trajectories can be 
found for which arrival occurs at Mercur)^ aphelion and 
when it is strongly inclined to the ecliptic. For (‘xample, 
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for a launch date of Septeinher 2, 1966, SPP 47.6 X 10 * 
km, LAP 6.39 deg and C, 109.9 krnVsec". It is 
interesting to note that tliesc* trajectory^ characteristics 
are almost duplicated for launch date August 26, 1973. 
Summarizing, it can he said that large variations in 
eeh'stial latitude and Sun-Mercury distance at encounter 
strongly influence' g(‘ocentric encagy, whc'reas the effc'ct 
of changes in helioccmtric central angle' is small. 

It should he noted that within the time interval eon- 
siderc'd in this report (1964-1975), Mercury aphelion 
occurs vc'ry nc'ar tlu' ecliptic so that trajectories with v('ry 
favorahic' c'lU'rgy characteristics can he found. As Mer- 
cury’s aphc'lion point moves away from the ecliptic, tlu' 
miniinum-encM'gy r('(|uir('ments incrc'ase. 

P"or Earth-Jiipiter transfers, the* geocentric energy will 
he smallest if (1) Jupiter is at perihelion (740 X 10' km) 
at encounter, (2) Jupiter is in the ecliptic at encountcT, 
(3) th(' heliocentric central angle is 180 deg, and (4) the 
Earth is at aphelion (147 X 10 ‘ km) at 1 aunch. 

Variations in absolute' minimum geocentric energy are 
caused hy tlie inclination (i ^ 1.3 deg) and eccentricity 
(e ^ 0.048) of jupih'r’s orhit and the elliptic orbit of 
the Earth. 

Tlu' traj('ctory with the' smalle'st energy in Tahk' 3-1 
(C ; 75.3 km Vs(‘C“) occurs for the launch date January 

3, 1970. Arrival occurs 985 days later, when Jupiter is 
close' te) tlu* ecliptic (LAP 0.004 deg) hut relative'ly far 
from perilielion (SPP 778.5 X KL km). 

Further study of Table* 3-1 shows that the strongest 
variations in geocentric e*ne'rgy are caused by large 
change*s in the Sun-Jupiter distance at arrival. In order to 
reduce* Earth-to-Jupiter flight times to two years, the 
geocentric energy reejuired is C, = 88 kmVsec" (launch 
date De'ce'mher 1, 1968). 

C, Mfssfon Payload 

The geocentric energy at launch C. determines the 
payload that can he placed on an interplanetary traje^c- 
tory. Figure 3-11 show^s approximate payload estimates 
vs Ci for three launch-\u*hicle configurations: two-stage 
Atlas/Centaur, three-stage Saturn and three-stage 
Saturn 5. These curvu's are only intended to give a gen- 
eral idea of payload; they were obtained from the open 
literature (Ref. 7). 

The maximum payload that can be placed on an 
Earth-Mercury trajectory using the Saturn IB launch 



Fig. 3-11. Typical payload capability curves for Atlas/ 
Centaur, Saturn IB, and Saturn 5 launch vehicles 


vehicle is 2500 lb. This payload is incr(*ased to 53,000 Ih 
if the Saturn 5 is used. 

Using the Saturn 5 vehich*, tlu* maximum payload that 
could he placed on an Earth-Jupiter trajectorv^ is 
26,000 Ih. 


D. Launch Period 

The practical design of a mission to Mercury and 
Jupiter will involvx* the establishment of a launch period, 
i.e., a period of consecutive days on which launch is 
possible. On each day the trajectories would have to 
satisfy numerous constraints, such as those placed on 
energy, communication distance, flight time, etc. For 
most missions, constraints on arrival conditions would 
exist also. 

As a first approximation, launch periods can be de- 
signed on the basis of geocentric energy, alone. Nat- 
urally, then, the minimum geocentric energy on ev^erv^ 
launch day is chosen. For each trajectory type the curves 
of minimum energy allow establishment of launch peri- 
ods of arbitrary length centered around each of the 
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absolute minimum-energy launch dates given in Table 
3-1. This selection has been carried out in Tables 3-2, 
3-3, and 3-4 for several of th(‘ Earth-Mercury and 
Earth-Jupiter opportunities tabulat(‘d in Table 3-1. 

For each trajectory type and class, launch periods of 
15, 30, and 45 days were determined. Eacli launch period 
has a maximum geocentric energy associated with it. 
This maximum value is reached on the first and last day 
of each launch period and results in the lowest payload. 

The minimum injection energy during the launch period 
can be found from the corresponding absolute minimum- 
energy trajectory in Table 3—1. It is of interest for orbiter 
missions since more fuel can be added to the retro tanks 
as the geocentric energy decreases during the launch 
period, resulting in closer orbits at the target planet. 
Table :3-2 also shows the variations, during the launch 
period, of the most important trajectory characteristics 
on which constraints are usually placed. 

The Earth- to-target-planet flight times Tr will affect 
virtually every spacecraft subsystem. The communication 
distance at arrival R( indicates the design of the com- 
munications system. The planetocentric asymptotic speed 
is a measure of the speed of the spacecraft at the 
time of closest approach to the target and is important 
for both orbiter and flyby missions. 

The geocentric asymptotic declination is included in 
Table 3-2 because it is an important parameter in de- 
termining the injection location over the Earth s surface. 
Acceptable values for this parameter are estimated to 
lie between —34 and +34 deg. Values outside this range 
result in severe restrictions for Cape Kennedy launchings. 

As an example, note how the 45-day launch period 
November 21, 1968 to January 5, 1969, for Type II, Class 
I trajectories to Jupiter was found from the curve of 
minimum geocentric energy. Fig. 3-10. Starting with the 
launch date December 13, 1968 (C;, — 77.8 kmVsec-), 
C:{ w'as increased until the range of launch dates was 45 
days. These limits then defined the launch period and 
the corresponding maximum energy (C.. ~ 86.2 kmVsec“) 
ov^er the period. 

Sometimes it happens that the curx^es of minimum 
geocentric energy corresponding to Type I and Type II 
trajectories intersect within the launch period for a given 
type. It then becomes possible to design launch periods 
during which both Type I and Type II trajectories would 
be used. The corresponding geocentric energies are lower 
than if either type were used alone. 


E. Generol Cbaracferisfks of Trajectories 

In order to gain further insight into the nature of 
Earth-Merc ury and Earth-Jupiter trajectories, plots of 
sevc'ral key trajectory variables were prepared for each 
of the launch opportunities. The relevant trajectory 
parameters can be categorized into three groups: geo- 
centric, heliocentric, and planetocentric. Sections IV-XV 
show detailed plots of the following trajectory param- 
c'ters, which are italicized: 

1. Geocentric Parameters 

During the geocentric phase, the magnitude and direc- 
tion of tlu^ hyperbolic-excess velocity vector V^p, or the 
velocity \a‘ctor of the spacecraft relative to the Earth at 
the time of injection into heliocentric orbit, are of prin- 
cipal interest. 

The magnitude of \hi> is represented by the geocentric 
energy at hunch, C , ^ All plots of trajectory param- 
eters show contours of constant and integral values of C j. 

The dir('( tion of also called the outgoing asymptote 
S, is given relative to several coordinate systems. The 
right ascen.sUm 0,s and declination <Ij,s of the outgoing 
asymptote specify its direction relative to a geocentric 
equatorial coordinate system. The celestial latitude y/, of 
S or the an^^/c between the geocentric asymptote and the 
launch phnefs orbital plane and the angle between the 
Sun-Eaiih vector and outgoing geocentric asymptote t/. 
specify the orientation of S relative to the heliocentric 
ecliptic coordinate system. 

Obvi(msly, the energy has a direct influence on space- 
craft weight. The declination of the outgoing asymptote 
also severely restricts the interplanetary trajectories which 
are possible in a more complicated way. 

2. Heliocentric Parameters 

The follow in g properties of heliocentric transfer ellipse 
are presented: aphelion distance Ra, perihelion distance 
Rp, and inclmation i with respect to the ecliptic. The 
eccentricity e and semimajor axis a of the heliocentric 
ellipse can be found from these. Sections IV-XV show 
graphs of true anomaly in the transfer ellipse at launch 
V i, and true anomaly in the transfer ellipse at arrival Vp, 
as well as of their difference, the heliocentric central 
angle +. 

The Earth-to-target-planet time of flight Tr and com- 
munication distance Rv are of direct interest to designers 
of all spacecraft subsystems. 
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Table 3-3. Type I Jupiter transfer characteristics 



4 / 5 - 4 / 20/73 15 0.8680 















Table 3-4. Type II Jupiter transfer characteristics 



5 / 4 - 5 / 19/73 15 0.8412 1329 1 475 630 
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The celestial latitude at arrival time I3p is referred to 
the ecliptic plane. 

3. Planetocentric Parameters 

The direction of the approach hyperbolic-excess velocity 
vector or velocity of the spacecraft with respect to 
the target planet at arrival, is described in several coor- 
dinate systems. 

The ri^ht ascension 0;, and declination of the unit 
vector along also known as the approach or incom- 
ing asymptote S, are referred to the Earth's ecpiatorial 
plane. The bettceen the incomin<^ asymptote and 

the arrival plaiiefs orbital plane and the angle bettceen 
the target planet-Sun vector and the incoming-planet- 
centered asymptote f;, relate the incoming asymptote 
direction to the orbital plane of the target planet. The 
angle between the planet-Earth vector and the incoming 
asymptote Ck defines the motion of the spacecraft relative 
to the Earth-planet line a few days prior to encounter. 
The angle between the planet-Canopus vector and the 
incoming asymptote Cr is of importance because Canopus 
will most likely be the second attitude reference for 
these missions, in addition to the Sun. 

The magnitude of the asymptotic speed with re- 
spect to the target planet is of great importance for 
orbiter, lander, and multiplanct missions. 

F. Discussion of Earfb-Mercury and 
Earth-Jupiter Trajectory Parameters^ 

1. Declination of the Geocentric Asymptote 

Equation (2-60) defined an iiuHpiality betwe<‘U 
(the launch azimuth), and <!>/, (the latitude of the 
launch site): 


The regions within which launch from Cape Kennedy 
can occur are shown in Fig. 3-12 as a function of 
Usually, the launch azimuth will he restricted to lie 
within certain bounds established by range safety con- 
siderations (typically 90 to 114 dc‘g). 

For the two Mercury opportunities considered (1967, 
1968), the declination constraint do(^s not restrict th(‘ 

•’A detailed discussion of the influence of the trajectory parameters 
on trajectory desij^n is pres<‘nted in Ret. 8 (S<‘ction Ill.D). The 
following material is intend(*d as snppletnentary to that prt'sented 
in Ref. 8 and 9. 



NORTH EAST SOUTH 

LAUNCH AZIMUTH 2^, deg 

Fig. 3—12. Permissible regions of the declination of the 
geocentric asymptote d^s for Cape Kennedy launchings 

Type I-(dass I trajcctoru\s, which have the shortest Right 
tiuK^s. Thv c'onstraint does eliminate a substantial part of 
the other trajj ctories (i.i\. Type I-Class 1 and Type 11- 
Class 1 and IF). This can be se(‘u from the corresponding 
plots of vs launch date (Figs. 4-5, 5-5). 

For tiu' 1968-69, 1969-70, and 1970-71 Jupiter oppor- 
tunities, th(‘ declination constraint does not restrict the 
possible T>’pc I -Class 1 trajectories (see Figs. 6-5, 7-5, 
8-5). In 1972 and 1973 it becomes almost impossible to 
list* this kiml of traj(*ctory (see Figs. 9—5, and 10-5). In 
fact, in 1972 and 1974 only the Type 11-Class 1 trajec- 
tories are not affected by the declination constraint. All 
other trajectory types are severely restricted by the 
constraint. 

2. Angle Between the Sun-Earth Vector and the 

Outgoing Geocentric Asymptote (i. 

This quantity, which is of importance for spacecraft 
using the Earth and Sun as attitude references, takes on 
values between 80 and 130 deg for the 1967 Mercury 
Type I trajectories (45 to 110 deg for Type II trajectories; 
see Fig. 4-8). It takes on approximately the same range 
of values during the 1968 Mercury opportunity owing 
to the fact that the Earth-Siin-Mercur>^ constellation is 
about the same for both opportunities. The quantity 
varies between 45 and 120 deg for the five Jupiter oppor- 
tunities considered herein (1968-69, 1969-70, 1970-71, 
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1972, 1973) for energies less than 120 kmVsec^ In all 
cases, it decreases with increasing launch date. 

3. Celestial Latitude of the Geocentric Asymptote y/. 

This quantity, referred to in the graphs as the angle 
between the outgoing geocentric asymptote and the 
launch planets orbital plane, is strongly dependent on 
the celestial latitude at arrival time (See Ref. 1, page 
35. Fig. 3-14 of Ref. 8 applies for Earth-Mercur\^ trajec- 
tories as well as for Earth-Venus trajectories.) 

For the Mercury opportunities in 1967 and 1968, the 
target planet is relatively close to the ecliptic at encoun- 
ter, which causes y/, to take on values on the order of 
10 deg. For Jupiter, y/, can take on values between ~45 
to -r45 deg for the energy range considered herein. Note 
that the Type I-Class I trajectories, which have the 
shortest flight times, all lie very close to the ecliptic. The 
Type Il-Class II trajectories, which correspond to the 
longest flight time, also show this characteristic. 

It might be asked how y/, can take on such high values 
in spite of the fact that Jupiter’s inclination is only 2.5 
deg. In Ref. 8, page 34, it is shown that 

V/, cos F/, sin i 
sin yi. = 77 

where V/, is the heliocentric velocity at launch, is the 
hyperbolic-excess velocity at launch, F/, is the heliocen- 
tric path angle at launch, and i is the inclination of the 
trajectory. The magnification factor V/7V/,^ is on the 
order of 4 for Earth-Jupiter trajectories. 

4. True Anomaly at Launch and Arrival and 

Heliocentric Transfer Angle 

In general, low-energy elliptic heliocentric orl)its from 
the Earth to the inner plan(4s (McTcairy, \T'iius) will 
have aphelion distance on the order of 1 an (astronomical 
unit) and perihelion distances on the order of the distance 
between the target planet at the Sun (0.3 an for Mer- 
cury). Conse(iuently, launch will occur near aph(4ion 
(true anomaly = 180 deg) and arrival near perihelion 
(true anomaly Vp — 0 deg). Figures 4-9 and 5-9 show 
that the variation in Vp is small, 160 to 200 deg for Type 
I and Type II Mercury trajectoric's. Note that launch 
occurs before aphelion at the beginning of the launch 
period, and after aphelion at the end of the launch period, 
as was the case for Earth-Vemis trajectories. 


Employing the same r(‘asoning, low-energy trajectories 
to the outer planets (Mars, Jupiter, etc.) liave pc‘rih(4ion 
distances on the order of 1 au and aphelion distanc(‘s on 
the order of the Sun-target planet distance. Figures 6-9, 

7- 9, 8-9, and 10-9 show that for an Earth-Jupiter tra- 
jectory, the true anomaly at launch remains within 30 
deg of perihelion (V/^ = 0 deg). Figures 6-10, 7-10, 

8- 10, 9-10, and 10-10 show that the Jupiter encountt*r 
can take place from about 40 deg b(dore to 30 d('g after 
aphelion. 

Figur(\s 4— through 10-3 show the variation in the 
heliocentric transfer angle or the difference betweiMi 
the true anomalies at launch and encounter, for Earth- 
Mercury and Earth-Jupiter trajectories, respectively. 
Note that none of the Type 1 or Type II contours inter- 
sect the “ 180 deg line. 

5. Aphelion and Perihelion of the Transfer Orbit 

Figures 4-12 and 5-12 show that for Earth-Mercury 
trajectories the aphelion distance (point farthest from the' 
Sun) varies from 1 au by less tlian 10%. The perihelion 
distances show stronger variations (Figs. 4-11 and 5-11). 
As mentioned in Ref. 8, the probe passes aphelion after 
launch if the true anomaly at launch is less than 180 deg; 
it passes perihelion before arrival (i.e., passes within 
Mercury’s orbit) if the true anomaly at arrival is positive. 

F’or Earth-Jupiter transfer*, the perihelion distance does 
not change drastically (Figs. 6- to 10-11). Note, how- 
ever, that the aphelion distances can become as large as 
2400 X 10'' km, or over 3 times the Sun-Jiipiter distance 
for geocentric energies C > ~ 120 km Vsec- (Figs. 6- 
througli 10-12). The plots of true anomaly show that 
aphelion woidd take place after Jupiter encounter. Owing 
to the exceptionally large mass of Jupiter, it will exert a 
very strong perturbation on the* trajeetory during en- 
counter. 

6, Time of Flight 

Typical Earth-Mercury flight times range from 80 to 
160 days, as is shown in Figs. 4-2 and .5-2. As the launch 
date progresses, it is seen that flight times are shortened 
for both Type I and Type II trajectories. Type II trajec- 
tories offer low('r energies and larger launch windows 
for a given (*nergy, with only a slight penalty in flight 
time. 
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Flight times to Jupiter are on tlie order of 1 to 4 years 
(Figs. C>- through 10-2). As launeh dat(’ is delayed, flight 
time increase's substantially for a gi\ en I'lic'rgy. 

7. Communication Distance 

For Mercury trajectorie's. tlu' communication distance 
at c'ncounter increase's as launch date' is delayeel, wliere'as 
flight time ek'crease's. The' e’ommunic.ition elistance* can 
take' on \'alue's be'twee'n 100 > 10' km to ISO X 10’' km 
(se'e Figs. 4—1 anel 5—1). 

The' communications distance's at jupite’r arrival e'xtend 
from 650 X 10’ km to 950 X 10" km. Since the Sun- 
Jupiter distance is 740 X 10" km and the Earth-Sun 


elistance' is 150 > 10" km, it is clenir that if the comrnuni- 

e'ations distance- at (‘iicounter is around 890 X 10" km, 
the' Sun will lin \ e-ry close* to the* Jupiter-Earth line' (see* 
k'igs. fX lhromi,h 10-4), 

8. Inclination of the Heliocentric Transfer Plane i and 

CX'lestial laitilude at Arrival 

The' inc lination i of the helioce'ntric transfe'r orbit to 
the* eeli[)tic jdane- is essentially a function of tw'o parain- 
e'tcTS, (1) the- Ik liocentric central angle T, and (2) the 
ee'le'Stial latitude* of the planc't at arri\al, since.' 


sin i 


sin /j;, 
sin ^ 



HELIOCENTRIC CENTRAL ANGLE, deg 

Fig. 3—13. Inclination of the heliocentric orbital plane to ecliptic vs heliocentric central angle 
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If tht‘ lu‘lioccntric central angle is fixc'd, and the absolute 
magnitudt' of the celestial latitude at arrival is increased, 
the inclination also increases. If the celestial latitude of 
the planet at arrival is zero, the inclination is equal to 
zero, no matter what the heliocentric ccaitral angle mav 
be (with the exception of a c(‘ntral angle of 180 deg). 
Howev(‘r, as shown in Fig. 3-13, if the celestial latitiuh* 

of the plaiK‘t at encounter is fixc'd at an absolute value 
greater than zero, and th(‘ ciaitral angle is varic'd, the 
inclination will reach its minimum valui' (equal to at 
central angh\s of 90 and 270 deg, or its maximum value of 
90 dc‘g at th(' central angles eciual to (3^,, 180 180 ^ 

and 360 -- l^i^. It is because of the fact that th(^ inclination 
increases near the central angles of 180 :hp,> that Type I- 
Class II and Type" Il-Class I trajectori(\s come into exist- 
ence. The incr(‘asing heliocentric inclination eventually 
results in a sharp increase in ener gy C , as th(" ccmtral angle" 
approaclu's 180 ele"g the‘reby bringing tlie* two groups 
of energy contours into existence. It may now become 
appare"ut tliat the inclination of the he*liocentric orbital 
plane* may take on values greater than the celestial lati- 
tude e)f the planet at arrival because of the varying helio- 
centric central angle. 

For the Me*rcury aiul Jupiter trajectories conside^red 
herein the inclinatie)n can be as high as 12 deg. By defi- 
nition, inclination is positive, whether encounter takes 
placi* above* or below the ecliptic. Figures 4-13 and 5-13 
show inclination vs launch date for Mercury. Figures 6- 
through 10-13 show this paramt'ter for Earth-Jupiter 
transfers. 

9. Asymptotic Approach Speed 

Figures 4-15 and 5-15 show that the variation in 
asymptotic speed with respect to Mercury is from 12 to 
24 km/sec for both Type I and Type II trajectories. Type 
I-Class I trajectories (shortest flight times) have higher 
approach speeds than Type I— Class II trajectories. It will 
be recalled that approach speeds to Venus were on the 
order of 3.75-10.7 km/sec. 

Approach speeds to Jupiter lie in the range of 5.5-14 
km/sec. The lowest asymptotic approach speeds are ob- 
tained using Type I-Class II or Type Il-Class I trajec- 
tories, shown in Figs. 6- through 10-15. 

10, Angle Between the Incoming Asymptote and the 

Arrival Planet s Orbital Plane yp 

For Mercury, Type I-Class I and Type Il-Class II 
trajectories have values of yp in the range —10 to +30 


deg. This angle can become as large as 60 deg for the 
other trajectories, as is explained in Ref. 8, page 51. (See 
Figs. 4-16 and 5-16.) For Jupiter trajectories, yp remains 
on the order of 15 deg or less for all opportunities. It 
takes on values on the order of 5 deg for Type I-Class I 
and Type Il-Class II trajectories. (See Figs. 6- through 
10-16.) 

11. Angle Between the Arrival Planet-Sun Vector and 
the Incoming Asymptote 

For Mercury Type 1 trajectories, the angle l^p is acute, 
whereas it varies between 60 and 160 deg for Type II 
trajectories. (See Figs. 4-17 and 5-17.) This angle is of 
considerable importance in the analysis of postencounter 
Jupiter trajectories (see Section III.F). As was shown for 
the Mars case in Ref. 8, page 51, f;, values for Type I- 
Class I trajectories are obtuse at the beginning of the 
launch windows and become acute at the end of the launch 
windows. The angle is smaller for the other Jupiter trajec- 
tories as is shown in Figs. 6 through 10-17. 

Figure 3-14 shows the illuminated crescent orientations 
at arrival for various values of Figures 3-15a and 3-15b 
show near-planet geometries for Type I and Type II 
approach trajectories to Mercury and Jupiter, respectively. 

Figures 3- 16a and 3-1 6b show these trajectories as 
viewed by an observer in inertial space. 

12. Right Ascension and Declination of the 
Incoming Planet-Centered Asymptote 

For Mercury trajectories, 0^, takes on values from 300 
to 355 deg. Right ascension 0^, ^ 360 deg corresponds to 
the vernal equinox direction, or direction of the Sun-Earth 
vector at the beginning of spring. The declination Op 
measured relative to the Earth's equator, varies from —40 
to 4 40 deg. 

For the Jupiter trajectories, there is a considerable 
variation in 0p for different launch opportunities: 160 to 
185 deg (1968-69), 290 to 320 deg (1973). The correspond- 
ing variations in 4>p are ~ 8 to 22 deg (1968-69) and —36 
to —8 deg (1973), respectively. 

13. Angle Between the Planet-Earth Vector and the 
Incoming Asymptote 

Plots of this parameter are not available for the Mercury 
trajectories at this time. The angle Ie is equivalent to the 
Earth-probe-target-planet angle a few days before arrival. 
The angle ~ 180 deg corresponds to motion along the 
Earth-target-planet line, away from the Earth, during the 
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(o)MERCURY TYPE I TRAJECTORY ( > 0 deg ) 

ENCOUNTER BEFORE PERIHELION 



(c^ MERCURY TYPE II TRAJECTORY ( > 0 deg) 

ENCOUNTER AFTER PERIHELION 




(b) MERCURY TYPE I TRAJECTORY ( y^ < Odeg) 
ENCOUNTER BEFORE PERIHELION 



(d) MERCURY TYPE n TRAJECTORY ( y^ < O deg ) 
ENCOUNTER AFTER PERIHELION 



<e) JUPITER TYPE I TRAJECTORY ( y^ > 0 deg ) 
ENCOUNTER BEFORE APHELION 


( f ) JUPITER TYPE I TRAJECTORY ( y^ < 0 deg) 
ENCOUNTER BEFORE APHELION 



(g) JUPITER TYPE I TRAJECTORY ( y^ > 0 deg) 
ENCOUNTER AFTER APHELION 


(h) JUPITER TYPE I TRAJECTORY ( y^ < 0 deg) 
ENCOUNTER AFTER APHELION 


Fig. 3—14. Crescent orientations for typical trajectories to Mercury and Jupiter 
as observed from approaching spacecraft 
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DIRECTION 
TO SUN 




(a) 



DIRECTION 
TO SUN 



(b) 


Fig. 3—15. Near-Jupiter geometry for typical trajectories, viewed from above ecliptic plane: (a) Type I trajectory, 
encountering Jupiter before aphelion; (b) Type I— II trajectory, encountering Jupiter after aphelion 


last (lays before' tar^ct-planct encounter; ~ 180 deg 
in('ans that tlu' prol>e will he oceiilte'd from view of th(^ 
Eartli ht'cause of hending of the trajectory hy the target 
planet. Figures 6- through 10-20 show that this param- 
eter vari(‘s hc'twt e'u 00 and 160 deg for Earth-jupiter 
trajectori(‘s. 

14. Angle Between the Planet-Canopus Vector and 
the Incoming Asymptote 

Plots of fr are not available for Earth-Mercury transfers 
at this time. For Earth-jupiter transfers, the ranges of fr 
are 75-150 deg (1968-69), 20-50 deg (1970-71), 0-20 deg 
(1972), 90-118 deg (1973), respectively; fr = 0 deg indi- 
catc's that the spacecraft will approach Jupiter along the 
spacecraft-Jupiter-Canopus line. (See Figs. 6- through 
10 - 21 .) 


G. Procedures for Utilization of Graphs in 
Design of Planetary Trajectories 

For a trajectory d('sign, a source of information must 
he availal)le which can he quickly scanned to determine 
tlie range of feasible trajectories for a mission. The graphs 
presemted in Sections IV to X of this report constitute 
such a source. The order and procedure for actual use of 
these graphs hy the trajectory engineer in the design and 
analysis of trajectories are now reviewed. 

(1) Both Type I and Type II transfers should be 
scanned. Those trajectories should be selected for 
which the declination of the outgoing geocentric 
asymptote lies roughly between —34 and {-34 deg. 
These are the feasible trajectories for launchings 
from AFETR. The algebraic value of the declina- 
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DIRECTION 
TO SUN 






Fig. 3—16. Near-Mercury geometry for typical trajectories, viewed from obove ecliptic plane: (a) Type I trajectory, 
encountering Mercury before perihelion; (bl Type II trajectory, encountering Mercury after perihelion 


tion rc'vcals vital inionnation concerning tlic' prc'- 
injc'ction trajc'etory (sc'c Section IlI.E.]), 

(2) For tlic trajc'ctorics of jxiragraph (1), the energy 
recpiireincaits (C ;) for various firing periods should 
l)e ol)servc'd. If it is found that booster payload 
capability and the' desired payload weight mateli 
the r(‘quired injection energy for a given firing 
peri()d (for example, 30 days) and also satisfy para- 
graph (1) for I'ype I and Type II trajectories, a 
dc^cision must tlum be made to i!tilize either Type I 
or Type II trajc^ctoric’s or, pcTliaps, both (sc'e 
Section III.D), 

(3) In making the decision to utili/e either the Type I 
or the Type II traj(‘ctory, or both, tlie curves of 
flight time and Earth-probe communication dis- 
tance vs launch date are most helpful. In general, 
Type II trajectoric‘s have longer flight times and 


Eai tli-pmbe distance's at (‘ucounte'r than do 1 > ])c I. 
I'he ae tual difh'rencc's in magnitude ch'j)end on the' 
mission and range of injc'ction (auagy. In geiK'ral, 
th(’ longc'r the* flight time' the greater is th(' seaisi- 
tivily of lhc‘ trajc'ctory te> injection c'n'ors. 

(4) \c‘xt. tlu' parameter f/., tlu' angle bclwaaMi th(’ 
outgoing as)in]:)tote and th<' Sun-Earth wclor at 
h(4io('( i itric inj('ction, should bc' stndi('(k It will b(' 
rec alled tliat this is cHjuivalc'nt to thc‘ I^arth-jirobe- 
Sun angle' at a fe'W days after launch. Since' tlu're 
aie usually many limitations on a sjXK'C’cralt whic h 
is stabili/c'd and controlled in attitude by optical 
reh'rc'nct's (such as the Eartlc Sun, and (ainopusb 
the' Isarth-]'>robe-Sun angle may bc' rc'stric'tc'd nesir 
the Earth aneb perhaps, throughout the' flight. TIm' 
parauK'ter is most lu'lpful in trajc'ctory design 
lor dc tninining the constraint 7U‘ar the Isarth. 
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(5) If tile Vtilu(‘ of f/, is less than 90 dcjf for Mercury 
traj('et()rit‘s, the prolx' is nsiially b(‘ing injected into 
tlu' helioecMitrie-transfer orbit Indore aph(9ion and 
will travc'l outside the Earth’s orbit before “falling 
in tow ard M{'rcury. To b(^ pR*cise, oiu* can observe 
tiu' cnrvf's of tru(' anomaly at lanneh. Rc‘cause of 
t('mp('ratiire'C‘ontrol probhans, tluax* may b(' a 
restriction on tlu' maximum Sun-probe distance dur- 
ing fliii;ht. I’o determim* this maximum distanc'C 
dnriinj; flight, om‘ simply finds tlu‘ aplu'Iion distance 
of tile prob(’ from the j^raphs. If is ^r(ait(a* than 
90 (h\u;, tlK‘ii tlie maximum distance^ durini^ ffi<j;ht to 
M(.a‘eur\ is {‘ss(Mitially tlu' Sun— prolu' distance at 
launth. 

(0) If the \ alu(‘ of is ^r(‘at(T than 90 d(*i^ for Jupiter 
trajc‘C“tori(‘s, th(‘ probe is usually bein^ injt‘cted into 
th(' lu9ioe(aitric-transfer orbit Indon* pca ihc'lion and 
w'ill “fall’’ insid(‘ tlu‘ Earth’s orbit Indore traversing 
out to )upit('r. TcanperatuR'-eontrol considcaations 
may rc'strict the minimum Sun-probe distanci^ dur- 
iiiL^ fliu;ht. To det(a-miiu‘ tlu’ minimum distanet\ 
sim]>ly find tlu‘ p(aih(‘lion of tin' piobc'’s transha* 
orbit for the di'siianl trajin tory oi' rant^e of trajen-- 
tori(vs. If Ci- than 90 d('^. tlu‘ minimum Sun- 

j)robc‘ distanc(‘ during; flight to [upit^a* is th(‘ 
distance' at launch. 

(7) For a sj:>ac“(‘C'raft stabilized and controlU'd in attitude* 
by Sun and Earth optical R*ft'r(*nc(*s, a constraint 
may c'xist which will R'strict the Earth-probc'-Sun 
angle* to a v alue* greaitea' than 0 de*g by a f(*w' degre'cs 
or ](*ss than ISO d(*g by a few' d(*gR*es. This imme*- 
diate'Iy implit*s that a trajec tory must be* c hose'u tliat 
has an inclination of the* he*lioee*ntri(‘ orbital j:)lane‘ 
to the* ('C'liptic' which is greater them 0 c/c'g. To satisfy 
the* above* constraint, the* inclination may have* to 
be 0.25, 0.50, or e*v(*n 1.0 de*g, depending on tlie* 
re'striction and on the* traje*ctory itsc'lf. The* valuevs 
of the* orbital inclinations for spe’C'ifie trajc'ctories 
can be* found from the* graphs in Se'ctions IV to X. 

(8) For c*sse*ntially all feasible I’ype* II trajc'ctoric's to 
Me*rc'ury and a few* of Type* I, Me*rcury (‘neounte*r 
for sorne* missions will take place affn' p(‘rihc*lion of 
the transfer orbit. This means that the probe will 
pass closest to the Sun s(*veral days or, perhaps, 
several wa*eks before encounter, dc'p(*nding on the 
trajectory. Be'cause of temperaturc'-control rc*(piire- 
rnents, there may be a R*striction on the minimum 
Sun-probe distance during flight. To determine* the 
value of this parame*ter, one finds the pc*rihelion 
distance of the probe for the desired trajectory from 


the* graphs. This paramc’tc'r is the* minimum Sun- 
])robe* distanc'c* during flight to Me*rc‘ury if the* tine 
anomaly of arrival is in the* first or se*cond epiadrant. 

(9) For Type* I traje*ctories to Jupiter, e*ncounte*r usually 
takc's ])lac(* be*fore* aphelion of the* transfer orbit. 
For those traje*ctorie*s having true anomaIie*s at 
arrival in the third or fourth (juadrant, however, 
c*ncounter will take* ]')lace* aftc*r aphe*lion. 1 his me*ans 
that the* maximum Sun-j)robe* distance during flight 
oc'curs be'fore* e*neonntc*r. d’empe*rature*-control con- 
side*ralions may re*stric*t the* maximum Sun-probe 
distance* during flight. To dc*te*rmine* this maximum 
for the* ne'cessarv^ tiaje*ctorie*s, find the* value* of the* 
aphe*lion distance from the* graphs in Se'ctions VI 
to X. 

(10) For plane'tary missions such as laude*rs or orbite*rs, 
it w ill probabK' be*come* ne*ee*ssary to utilize* traje'c- 
torie's tor which the* hyj')e*rbolic‘-e*xc'e*ss sj)e*e*d at the 
plane't is ne*arly minimum. In fact, in orde*r to 
maximize* the* scic*ntific-j:)ayload we*ight, a “tradeoff” 
must be* made* in minimizing the ge*oc*e*ntric e*ne'rg\' 
or planetoce*ntric e*ne*rgy. Suc'h traje'ctorie'S will 
allow the* he*a\'ie*st j')ayloads to be* lande*d on the* 
planet or inje*cte*d in a de*sire*d e*lliptical orbit around 
the* plane't w ith the* use of re*tro-mane'uve*rs. 9 he* 
traje'ctorv de*signe*r can find the* values of the* 
h\'[:>e*rbolic-e*xce’ss spe*e‘cls at the* plane*! for various 
traje‘C‘toric*s from the* graphs in Sc'ctions IV to X, and 
c*an dc*te*rmine* the* fe*asibility of the* traje*ctorie*s from 
the* ])i'oct*dure*s in ]:>aragraphs (1) to (9). 

(11) For a si)ac*e*c*iaf t which is stabiliz.c*cl and C'ontrolIe*d 
in attitude* by 0 ]')tical re*ference*s, the*re* may be* 
cc'iiain re'strictions on the Sun-probe-plane*t angle 
as the* probe* approaclu's the plane*!. Such re*stric'tions 
mav' ])e* ne*ce*ssary bc'cause* of a]')proach-guidance* 
conside'rations. The* ])aram<*te*r the* angle* bc*- 
tvve*(*n the* ai)]')roach asvinptote* and the targc't- 
j:)lane*t-Sun ve'c tor, is e*ejuivale*nt to the* Snn-]’)robe'- 
plane't angle* a fe*w days be*lore* e*ne'Ounte*r. The* 
l)arame*te*r csui be* most he*lpful in analyzing pe*r- 
missible* traje*ctorie*s ne*ar the* plane*!. 

(12) lU*garding the de'sign of traje*ctorie*s ne*ar the planet, 
it is n(*e*dle*ss to say that, by altering tlu* parameters 
in the* pi-e*inje*ction traje*ctory nc*ar Earth by small 
incre'inents, the* probe can be made* to pass on any 
side* of the* plane*t. Hovve*ver, the* inclinations of the 
plane'toc'entric hyperbolic orbit (or e*l!iptical orbit, 
assuming that the* re*(piire*d retro-maneuv e*r is made 


45 



JPL TECHNICAL REPORT NO. 32-77 


in tlu' orbital plaiu' ol the In prrbola) to tlie planet s 
orbital plant' dt'pt'iul.s on tlu' aiming point at the 
plant'b as wt'll as on the* param(‘t(*r y,.,, tlu‘ anglt* 
bet\\(*en the approadi asymptote atul the planets 
orbital plant'. The minimum inelination which ean 
bt' attained is t*(pial to 7^,. 11 y^, is zero, the inclina- 
tion of the orbit to tbt' planet’s orbital plant' will 
range' from 0 to ISO deg, depending on the aiming 
pt)int. Inelinations ol 0 to 90 dt'g irnply direct 
motion (in the dirt'clion ol the [:)Ianet s orbital rota- 
tion), \\ ht'reas inc linations from 90 to ISO deg imply 
i c'trogradc' motion (opposite' to tlu' ]')lanet s orbital 
i-otation). If y,, ‘ 45 deg, tlu' inelination may 

range' from 45 to 135 dc'g lor all aiming points. 
II yi, ' 90 tlt'g, the' inclination will t‘(|nal 90 dt'g 
lor anv aiming point ^^bic‘h is ehosi'ii at the planet. 
lOr 7,. ; > 0 dt'g, th(' [’)mbt' w ill approach the' plant't 
Irom a [)ath below' tht' i^lanct's orbital plant'; lor 
7^, < 0 dt'g, the' aj^])roach w ill bt' Irom above. In 
main' cases, it max' bt' dt'sirablc' to dt'sign the nt'ar- 
])laiu't trajectory with a prt'scribt'cl inclination in 
mind such that the jirobt' will [xiss the target body 
in the' plant't s orbital or, ]M'rhaj)s, (‘(piatorial plane. 
It is thus apj'iart'nt that the paramc'tc'r y^, max' or 
may not [H'rmit the' sc'lt'cted pass. To find the valnt' 
ol 7^, lor a gixt'n trajc'ctorx', obsc'rxc' the cnrx t' of y„ 
xs launch datt' lor tht' dt'sirt'd mission. 

(13) To comprt'hend Inlly tht' significance of t'ach 
parainc'tt'r and its variation xxitli launch date and 
arrix al date', it max' bt' adx antagt'ons to coirstrnct 
tht' loci ol constant arrixal dates on t'ach pertinent 
graph in Sectit)ns I\' to X. I his may bt' dt>nt' in 
three' stt‘[)s, as shoxx n in Fig. 3-17. 

(a) C'onstrnct tht' loci of dt'sirt'd arrix al dates t)n the 
curves of flight time' xs launch date, 

(1)) Mark the latmch date's and gc'ocentric energit's 
lor xvhich interst'ction ol tlu' arrixal loci and the 
xarions t'lu'rgy contours occur. 


STEPS 
I AND 2 



LAUNCH DATE 



LAUNCH DATE 


Fig. 3—17. Steps for construction of 
arrival date loci 


(c) Tak<.' th(' intt'rsection points of (b) and construct 
tht' ;n rixabdatc' loci on the graph ol intc'rest. 


Ex])t'i i('nc'( has shoxvn that trajc'ctory parameters tend 
to t'xiabit an inxariance if tlu' arrival date is ht*ld fixed 
XX hilt* tin' launch date is xarit'd. This fact has provt'd vt'ry 
nst'lul in trajt'C'tory design. 
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IV. MERCURY 1967: TRAJECTORY PARAMETER GRAPHS 

Figure 

4-1. Mercury 1967: Minimum injection energy vs launch date 

4—2. Mercury 1967: Time of flight vs launch date 

4-3. Mercury 1967: Heliocentric central angle vs launch date 

4—4. Mercury 1967: Earth-Mercury communication distance vs 
launch date 

4—5(1). Mercury 1967: Declination of the geocentric asymptote vs 
launch date, Type I 

4—5(11). Mercury 1967: Declination of the geocentric asymptote vs 
launch date, Type II 

4—6(1). Mercury 1967: Right ascension of the geocentric asymptote vs 
launch date, Type I 

4—6(11). Mercury 1967: Right ascension of the geocentric asymptote vs 
launch date, Type II 

4—7(1). Mercury 1967: Angle between outgoing geocentric asymptote 
and launch planet’s orbital plane vs launch date, Type I 

4—7(11), Mercury 1967: Angle between outgoing geocentric asymptote 
and launch planet's orbital plane vs launch date, Type II 

4—8(1). Mercury 1967: Angle between Sun— Earth vector and outgoing 
geocentric asymptote vs launch date, Type I 

4—8(11). Mercury 1967: Angle between Sun— Earth vector and outgoing 
geocentric asymptote vs launch date. Type II 

4—9(1). Mercury 1967; True anomaly in transfer ellipse at launch 
time vs launch date. Type I 

4-9(11). Mercury 1967; True anomaly in transfer ellipse at launch 
time vs launch date, Type II 

4—10. Mercury 1967: True anomaly in transfer ellipse at arrival 
time vs launch date 

4—1 1(1). Mercury 1967: Perihelion of transfer orbit vs launch date, Type I 

4—11(11). Mercury 1967: Perihelion of transfer orbit vs launch date, Type II 

4—12(1). Mercury 1967: Aphelion of transfer orbit vs launch date. Type I 

4—12(11). Mercury 1967: Aphelion of transfer orbit vs launch date. Type II 

4—13(1). Mercury 1967: Inclination of the heliocentric transfer plane 
vs launch date, Type I 

4—13(11). Mercury 1967: Inclination of the heliocentric transfer plane 
vs launch date. Type II 

4—14. Mercury 1967: Celestial latitude at arrival time vs launch date 
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IV. MERCURY 1967: TRAJECTORY PARAMETER GRAPHS (Cont’d) 


Figure 

4-15(1). Mercury 1967: Asymptotic speed with respect to Mercury vs 
launch date, Type I 

4-15(11). Mercury 1967: Asymptotic speed with respect to Mercury vs 
launch date, Type II 

4—16(1). Mercury 1967: Angle between incoming hermiocentric 
asymptote and arrival planet’s orbital plane vs 
launch date, Type I 

4—16(11). Mercury 1967: Angle between incoming hermiocentric 
asymptote and arrival planet’s orbital plane vs 
launch date, Type II 

4_17, Mercury 1967: Angle between Mercury-Sun vector and incoming 
hermiocentric asymptote vs launch date 

4-18(1). Mercury 1967: Declination of the hermiocentric asymptote vs 
launch date, Type I 

4-18(11). Mercury 1967: Declination of the hermiocentric asymptote vs 
launch date. Type II 

4-19(1). Mercury 1967: Right ascension of hermiocentric asymptote vs 
launch date, Type I 

4-19(11). Mercury 1967: Right ascension of hermiocentric asymptote vs 
launch date, Type II 
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f'9- 4—1. Mercury 1967: Minimum injection energy vs launch date 
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Fig. 4-3. Mercury 1967: Heliocentric central angle vs launch date 
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Mercury 1967: Earth— Mercury communication distance vs launch date 
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LAUNCH DATE 

Fig. 4—5(1). Mercury 1967: Declination of the geocentric asymptote vs launch date, Type I 
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LAUNCH DATE 

Fig. 4-5(11). Mercury 1967: Declination of the geocentric asymptote vs launch date. Type II 
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Fig. 4—6(1). Mercury 1967: Right ascension of the geocentric asymptote vs launch date, Type I 
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Fig. 4—7(1). Mercury 1967: Angle between outgoing geocentric asymptote and launch planet's orbital plane vs launch date, Type I 
















Fig. 4—8(1). Mercury 1967: Angle between Sun— Earth vector and outgoing geocentric asymptote vs launch date^ Type 1 
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Fig. 4—9(1). Mercury 1967: True anomaly in transfer ellipse at launch time vs launch date, Type I 
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10. Mercury 1967: True anomaly in transfer ellipse at arrival time vs launch date 
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1 1(1). Mercury 1967: Perihelion of transfer orbit vs launch date, Type I 
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Fig. 4—1 1(11). Mercury 1967: Perihelion of transfer orbit vs launch date, Type II 
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1—12(1). Mercury 1967; Aphelion of transfer orbit vs launch date, Type I 
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Fig. 4-13(1). Mercury 1967: Inclination of the heliocentric transfer plane vs launch date, Type I 
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Fig. 4-13(11). Mercury 1967: Inclination of the heliocentric transfer plane vs launch date, Type II 
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14. Mercury 1967: Celestial latitude at arrival time vs launch date 
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Fig. 4—15(1). Mercury 1967: Asymptotic speed with respect to Mercury vs launch date, Type I 
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LAUNCH DATE 

16(11). Mercury 1967: Angle between incoming hermiocentric asymptote and arrival planet’s orbital plane vs launch date, Type II 










Fig. 4—17. Mercury 1967: Angle between Mercury— Sun vector and incoming hermiocentric asymptote vs launch date 




























Pig, 4-18(1). Mercury 1967: Declination of the hermiocentric asymptote vs launch date, Type I 






























CLASS I 
CLASS n 





























LAUNCH DATE 

19(1). Mercury 1967: Right ascension of hermiocentric asymptote vs launch date, Type I 
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V. MERCURY 1968: TRAJECTORY PARAMETER GRAPHS 

Figure 

5-1. Mercury 1968: Minimum injection energy vs launch date 

5—2. Mercury 1968: Time of flight vs launch date 

5-3. Mercury 1968: Heliocentric central angle vs launch date 

5_4, Mercury 1968: Earth-Mercury communication distance vs 
launch date 

5-5(1). Mercury 1968: Declination of the geocentric asymptote vs 
launch date, Type I 

5-5(11), Mercury 1968: Declination of the geocentric asymptote vs 
launch date, Type II 

5-6(1). Mercury 1968: Right ascension of the geocentric asymptote 
vs launch date, Type I 

5-6(11). Mercury 1968: Right ascension of the geocentric asymptote 
vs launch date, Type II 

5_7(|). Mercury 1968: Angle between outgoing geocentric asymptote 
and launch planet’s orbital plane vs launch date, Type I 

5-7(11). Mercury 1968: Angle between outgoing geocentric asymptote 
and launch planet’s orbital plane vs launch date, Type II 

5-8(1). Mercury 1968: Angle between Sun-Earth vector and outgoing 
geocentric asymptote vs launch date. Type I 

5-8(11). Mercury 1968: Angle between Sun-Earth vector and outgoing 
geocentric asymptote vs launch date, Type II 

5_9(|). Mercury 1968: True anomaly in transfer ellipse at launch time 
vs launch date, Type I 

5-9(11). Mercury 1968: True anomaly in transfer ellipse at launch time 

vs launch date. Type II 

5-10. Mercury 1968: True anomaly in transfer ellipse at 
launch time vs launch date 

5-11(1). Mercury 1968: Perihelion of transfer orbit vs 
launch date. Type I 

5-1 1(11). Mercury 1968: Perihelion of transfer orbit vs 
launch date, Type II 

5-12(1). Mercury 1968: Aphelion of transfer orbit vs 
launch date, Type I 

5-12(11), Mercury 1968: Aphelion of transfer orbit vs 
launch date. Type II 

5-13(1). Mercury 1968: Inclination of the heliocentric transfer plane 
vs launch date, Type I 
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V. MERCURY 1968: TRAJECTORY PARAMETER GRAPHS (Cont’d) 


Figure 

5—13(11). Mercury 1968: Inclination of the heliocentric transfer plane 
vs launch date, Type II 

5—14. Mercury 1968: Celestial latitude at arrival time vs launch date 

5-15(1). Mercury 1968: Asymptotic speed with respect to Mercury vs 
launch date, Type I 

5-15(11). Mercury 1968: Asymptotic speed with respect to Mercury vs 
launch date, Type II 

5—16(1). Mercury 1968: Angle between the incoming hermiocentric 
asymptote and arrival planet's orbital plane vs 
launch date, Type I 

5—16(11). Mercury 1968: Angle between the incoming hermiocentric 
asymptote and arrival planet's orbital plane vs 
launch date, Type II 

5—17. Mercury 1968: Angle between Mercury— Sun vector and incoming 
hermiocentric asymptote vs launch date. Type I 

5-18(1). Mercury 1968: Declination of the hermiocentric asymptote vs 
launch date. Type I 

5-18(11). Mercury 1968: Declination of the hermiocentric asymptote vs 
launch date, Type II 

5—19(1). Mercury 1968: Right ascension of the hermiocentric asymptote 
vs launch date, Type I 

5-19(11). Mercury 1968: Right ascension of the hermiocentric asymptote 
vs launch date, Type II 
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Fig. 5-3. Mercury 1968: Heliocentric central angle vs launch date 




























JPL TECHNICAL REPORT NO. 32-77 



69p ‘*( 1 ) aiOidWASV 0iaiN30039 3H1 30 N0llVNn03a 


Fig. 5—5(1). Mercury 1968: Declination of the geocentric asymptote vs launch date, Type I 
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Fig. 5—5(11). Mercury 1968: Decimation of the geocentric asymptote vs launch date^ Type II 


























LAUNCH DATE 

Fig. 5-6(1). Mercury 1968: Right ascension of the geocentric asymptote vs launch date, Type I 































Fig. 5—6(11). Mercury 1968: Right ascension of the geocentric asymptote vs launch date, Type II 
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Fig. 5-7(1). Mercury 1968: Angle between outgoing geocentric asymptote and launch planet’s orbital plane vs launch date, Type I 




Fig. 5—7(11). Mercury 1968: Angie between outgoing geocentric asymptote and launch planet’s orbital plane vs launch date, Type II 









































Fig. 5-8(1). Mercury 1968: Angle between Sun-Earth vector and outgoing geocentric asymptote vs launch date, Type I 
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Fig. 5—8(11). Mercury 1968: Angle between Sun— Earth vector and outgoing geocentric asymptote vs launch date, Type II 
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Fig. 5—12(1). Mercury 1968: Aphelion of transfer orbit vs launch date, Type I 
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Fig. 5—13(1). Mercury 1968: Inclination of the heliocentric transfer plane vs launch date, Type I 













Fig. 5-13(11). Mercury 1968: Inclination of the heliocentric transfer plane vs launch date, Type II 
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14. Mercury 1968: Celestial latitude at arrival time vs launch date 
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15(l). Mercury 1968: Asymptotic speed with respect to Mercury vs launch date, Type I 
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15(11). Mercury 1968: Asymptotic speed with respect to Mercury vs launch date, Type II 



















16(11. Mercury 1968: Angle between the incoming hermiocentric asymptote and arrival planet’s orbital plane vs launch date, Type I 
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Fig. 5“16(ll). Mercury 1968: Angle between the incoming hermiocentric asymptote and arrival planet’s orbital plane vs launch date, Type II 
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17. Mercury 1968: Angle between Mercury— Sun vector and incoming hermiocentric asymptote vs launch date, Type I 
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Fig. 5—18(1). Mercury 1968: Declination of the hermiocentric asymptote vs launch date, Type I 
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18(11). Mercury 1968: Declination of the hermiocentric asymptote vs launch date, Type il 
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VL JUPITER 1968-69: TRAJECTORY PARAMETER GRAPHS 


Figure 

6-1. 

6-2. 

6-3. 

6-4(1). 

6-4(11). 

6-5(1). 

6-5(11). 

6 - 6 ( 1 ). 

6 - 6 ( 11 ). 

6-7(1). 

6-7(11). 

6 - 8 ( 11 . 

6 - 8 ( 11 ). 

6-9(1). 

6-9(11). 

6 - 10 . 

6 - 11 ( 1 ). 

6 - 11 ( 11 ) 

6 - 12 ( 1 ). 


Jupiter 68-69: Minimum injection energy vs launch cJate 

Jupiter 68—69: Time of flight vs lounch date 

Jupiter 68—69: Heliocentric central angle vs launch date 

Jupiter 68—69: Earth-Jupiter communication distance vs 
launch date. Type I 

Jupiter 68-69: Earth-Jupiter communication distance vs 
launch date, Type II 

Jupiter 68—69: Declination of the geocentric asymptote vs 
launch date, Type I 

Jupiter 68—69: Declination of the geocentric asymptote vs 
launch date. Type II 

Jupiter 68—69: Right ascension of the geocentric osymptote 
vs launch date, Type I 

Jupiter 68-69: Right ascension of the geocentric asymptote 
vs launch date. Type II 

Jupiter 68—69: Angle between the outgoing geocentric 
asymptote and launch planet's orbital plane vs 
launch date. Type I 

Jupiter 68-69: Angle between the outgoing geocentric 
asymptote and launch planet’s orbital plane vs 
launch date, Type II 

Jupiter 68—69: Angle between Sun— Earth vector and the 
outgoing geocentric asymptote vs launch date. Type I 

Jupiter 68-69: Angle between Sun— Earth vector and the 
outgoing geocentric asymptote vs launch date, Type II 

Jupiter 68—69: True anomaly In transfer ellipse at 
launch time vs launch date. Type I 

Jupiter 68-69: True anomaly in transfer ellipse at 
launch time vs launch date. Type II 

Jupiter 68—69: True anomaly in transfer ellipse at 
arrival time vs launch date 

Jupiter 68—69: Perihelion of transfer orbit vs launch date. 
Type I 

Jupiter 68-69: Perihelion of transfer orbit vs launch date. 

Type II 

Jupiter 68—69: Aphelion of transfer orbit vs launch date, 
Type I 
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VI. JUPITER 1968-69: TRAJECTORY PARAMETER GRAPHS (Cont’d) 


figure 

6 - 12 ( 11 ). 

6-13(1). 

6-13(11). 

6-14. 

6-15(1). 

6-15(11). 

6-16. 

6-17. 

6-18(1). 

6-18(11). 

6-19(1). 

6-19(11). 

6 - 20 ( 1 ). 

6 - 20 ( 11 ). 

6 - 21 ( 1 ). 

6 - 21 ( 11 ). 


Jupiter 68—69: Aphelion of transfer orbit vs launch date, 

Type II 

Jupiter 68—69: Inclination of the heliocentric transfer 
plane vs launch date, Type I 

Jupiter 68—69: Inclination of the heliocentric transfer 
plane vs launch date, Type II 

Jupiter 68—69: Celestial latitude at arrival time vs launch date 

Jupiter 68—69: Asymptotic speed with respect to Jupiter vs 
launch date, Type I 

Jupiter 68—69: Asymptotic speed with respect to Jupiter vs 
launch date, Type II 

Jupiter 68—69: Angle between incoming zeocentric asymptote 
and arrival planet’s orbital plane vs launch date 

Jupiter 68—69: Angle between Jupiter— Sun vector and incoming 
zeocentric asymptote vs launch date 

Jupiter 68—69: Declination of the zeocentric asymptote vs 
launch date, Type I 

Jupiter 68—69: Declination of the zeocentric asymptote vs 
launch date. Type II 

Jupiter 68—69: Right ascension of zeocentric asymptote vs 
launch date, Type I 

Jupiter 68-69: Right ascension of zeocentric asymptote vs 
launch date, Type II 

Jupiter 68—69: Angle between Planet— Earth vector and 
incoming asymptote vs launch date, Type I 

Jupiter 68—69: Angle between Planet— Earth vector and 
incoming asymptote vs launch date, Type II 

Jupiter 68-69: Angle between the planet-Canopus vector 
and incoming asymptote vs launch date, Type I 

Jupiter 68—69: Angle between the planet-Canopus vector 
and incoming asymptote vs launch date, Type II 
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Fig. 6—3. Jupiter 68—69: Heliocentric central angle vs launch date 
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Fig. 6—4(1). Jupiter 68—69: Earth— Jupiter communication distance vs launch date, Type I 






























0001 


JPL TECHNICAL REPORT NO. 32-77 



119 


Fig. 6— 4{ll). Jupiter 68—69: Earth— Jupiter communication distance vs launch date^ Type II 
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Fig. 6—5(11). Jupiter 68—69: Declination of the geocentric asymptote vs launch date, Type II 
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Fig. 6—6(1). Jupiter 68—69: Right ascension of the geocentric asymptote vs launch (date, Type I 



220 


JPL TECHNICAL REPORT NO, 32-77 



123 


Fig. 6—6(11). Jupiter 68—69: Right ascension of the geocentric asymptote vs launch date, Type II 












7(1). Jupiter 68—69: Angle between the outgoing geocentric asymptote and launch planet's orbital plane vs launch date, Type I 
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Fig. 6— 8(1). Jupiter 68—69: Angle between Sun— Earth vector and the outgoing geocentric asymptote vs launch date, Type I 
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Fig. 6—8(11). Jupiter 68—69: Angle betvy/een Sun— Earth vector and the outgoing geocentric asymptote vs launch date, Type II 
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Fig. 6-10. Jupiter 68-69: True anomoly in transfer ellipse at arrival time vs launch date 
















1 1(1). Jupiter 68—69: Perihelion of transfer orbit vs launch date, Type I 
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Fig. 6—1 1 (ID. Jupiter 68—69: Perihelion of transfer orbit vs launch date, Type II 
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12(11). Jupiter 68—69: Aphelion of transfer orbit vs launch date, Type II 
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14. Jupiter 68—69: Celestial latitude at arrival time vs launch date 









15(1). Jupiter 68—69: Asymptotic speed with respect to Jupiter vs launch date, Type I 
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16. Jupiter 68—69: Angle between incoming zeocentric asymptote and arrival planet’s orbital plane vs launch date 
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17. Jupiter 68—69: Angle between Jupiter— Sun vector and incoming zeocentric asymptote vs launch date 
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18(1). Jupiter 68—69: Declination of the zeocentric asymptote vs launch date, Type I 
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Fig. 6—18(11). Jupiter 68—69: Declination of the zeocentric asymptote vs launch date, Type II 
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Fig. 6—19(1). Jupiter 68—69: Right ascension of zeocentric asymptote vs launch date, Type I 
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1-20(1). Jupiter 68-69: Angle between Planet-Earth vector and incoming asymptote vs launch date, Type I 























Fig. 6—20(11). Jupiter 68—69: Angle between Planet— Earth vector and incoming asymptote vs launch date, Type II 
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Fig. 6-21(1). Jupiter 68-69: Angle between the planet-Canopus vector and incoming asymptote vs launch date, Type I 
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1—21(11). Jupiter 68—69: Angle between the planet— Canopus vector and incoming asymptote vs launch date, Type II 
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VII. JUPITER 1969-70: TRAJECTORY PARAMETER GRAPHS 

Figure 

7—1. Jupiter 69-70: Minimum injection energy vs launch date 

7-2. Jupiter 69-70: Time of flight vs launch date 

7-3. Jupiter 69-70: Heliocentric central angle vs launch date 

7—4(1). Jupiter 69—70: Earth— Jupiter communication distance vs 
launch date, Type I 

7-4(11). Jupiter 69-70: Earth-Jupiter communication distance vs 
launch date, Type II 

7-5(1). Jupiter 69-70: Declination of the geocentric asymptote vs 
launch date, Type I 

7—5(11). Jupiter 69-70: Declination of the geocentric asymptote vs 
launch date, Type II 

7-6(1). Jupiter 69-70: Right ascension of the geocentric asymptote vs 
launch date, Type I 

7—6(11). Jupiter 69-70: Right ascension of the geocentric asymptote vs 
launch date, Type II 

7-7(1). Jupiter 69-70: Angle between outgoing geocentric asymptote 
and launch planet's orbital plane vs launch date, Type I 

7-7(11). Jupiter 69-70: Angle between outgoing geocentric asymptote 
and launch planet's orbital plane vs launch date, Type II 

7—8. Jupiter 69— 70: Angle between Sun— Earth vector and outgoing 
geocentric asymptote vs launch date 

7—9. Jupiter 69—70: True anomaly in transfer ellipse at launch 
time vs launch date 

7-10. Jupiter 69-70: True anomaly in transfer ellipse at arrival 
time vs launch date 

7—11(1). Jupiter 69— 70: Perihelion of transfer orbit vs launch date. 

Type I 

7-11(11). Jupiter 69-70: Perihelion of transfer orbit vs launch date, 

Type II 

7—12(1). Jupiter 69— 70: Aphelion of transfer orbit vs launch date. 

Type I 

7—12(11). Jupiter 69— 70: Aphelion of transfer orbit vs launch date, 

Type II 

7-13(1). Jupiter 69-70: Inclination of the heliocentric transfer plane 
vs launch date, Type I 

7-13(11). Jupiter 69-70; Inclination of the heliocentric transfer plane 
vs launch date. Type II 
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VII. JUPITER 1969-70: TRAJECTORY PARAMETER GRAPHS (Cont’d) 


Figure 

7-14. 

7-15(1). 

7-15(11). 

7-16(1). 

7-16(11). 

7-17. 

7-18(1). 

7-18(11). 

7-19(1). 

7-19(11). 


Jupiter 69—70: Celestial latitude at arrival time vs launch date 

Jupiter 69—70: Asymptotic speed with respect to Jupiter 
vs launch date, Type 1 

Jupiter 69—70: Asymptotic speed with respect to Jupiter 
vs launch date, Type II 

Jupiter 69—70: Angle between incoming zeocentric asymptote 
and arrival planet's orbital plane vs launch date, Type 1 

Jupiter 69—70: Angle between incoming zeocentric asymptote 
and arrival planet's orbital plane vs launch date, Type II 

Jupiter 69-70: Angle between Jupiter-Sun vector and 
incoming zeocentric asymptote vs launch date 

Jupiter 69—70: Declination of the zeocentric asymptote vs 
launch date, Type I 

Jupiter 69—70: Declination of the zeocentric asymptote vs 
launch date, Type II 

Jupiter 69—70: Right ascension of the zeocentric asymptote vs 
launch date, Type I 

Jupiter 69—70: Right ascension of the zeocentric asymptote vs 
launch date. Type II 
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Fig. 7—3. Jupiter 69—70: Heliocentric central angle vs launch date 
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Fig, 7—4(1). Jupiter 69-70: Earth-Jupiter communication distance vs launch date, Type I 
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Fig. 7—5(1). Jupiter 69—70: Declination of the geocentric asymptote vs launch date, Type I 
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Fig. 7—5(11). Jupiter 69—70: Declination of the geocentric asymptote vs launch date, Type II 
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Fig. 7—6(1). Jupiter 69—70: Right ascension of the geocentric asymptote vs launch cfate, Type I 
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Fig. 7-7(1). Jupiter 69-70: Angle between outgoing geocentric asymptote and launch planet’s orbital plane vs launch date, Type I 
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-7(11). Jupiter 69-70: Angle between outgoing geocentric asymptote and launch planet’s orbital plane vs launch date, Type II 
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Fig. 7—8. Jupiter 69—70: Angle between Sun— Earth vector and outgoing geocentric asymptote vs launch date 
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Fig. 7—11(11). Jupiter 69—70: Perihelion of transfer orbit vs launch date, Type It 
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Fig. 7—12(11). Jupiter 69—70: Aphelion of transfer orbit vs launch date. Type II 
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13(11). Jupiter 69—70: Inclination of the heliocentric transfer plane vs launch date. Type II 


































14. Jupiter 69—70: Celestial latitude at arrival time vs launch date 
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18(1). Jupiter 69-70; Declination of the zeocentric asymptote vs launch date, Type I 
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— 18(11). Jupiter 69—70; Declination of the zeocentric asymptote vs launch date, Type II 
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19(1). Jupiter 69-70; Right ascension of the zeocentric asymptote vs launch date, Type I 
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19(11). Jupiter 69—70: Right ascension of the zeocentric asymptote vs launch date, Type II 
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VIII. JUPITER 1970-71: TRAJECTORY PARAMETER GRAPHS 

Figure 

8—1. Jupiter 70— 71 : Minimum injection energy vs launch date 

8-2. Jupiter 70-71 : Time of flight vs launch date 

8—3. Jupiter 70— 71 : Heliocentric central angle vs launch date 

8-4(1). Jupiter 70— 7 1 : Earth— Jupiter communication distance vs 

launch date. Type I 

8— 4111). Jupiter 70— 71 : Earth— Jupiter communication distance vs 
launch date, Type il 

8—5(1). Jupiter 70— 71 : Declination of the geocentric asymptote vs 

launch date, Type I 

8-5(11). Jupiter 70-71 : Declination of the geocentric asymptote vs 
launch date, Type II 

8—6. Jupiter 70— 71 : Right ascension of the geocentric asymptote vs 

launch date 

8—7. Jupiter 70— 71 : Angie between outgoing geocentric asymptote 

and launch planet's orbital plane vs launch date 

8—8(1). Jupiter 70— 71 : Angle between Sun— Earth vector and outgoing 
geocentric asymptote vs launch date, Type I 

8-8(11). Jupiter 70— 71 : Angle between Sun— Earth vector and outgoing 
geocentric asymptote vs launch date. Type II 

8—9(1). Jupiter 70— 71 : True anomaly in transfer ellipse at launch 
time vs launch date, Type I 

8—9(11). Jupiter 70— 71 : True anomaly in transfer ellipse at launch 
time vs launch date, Type II 

8-10. Jupiter 70— 71 : True anomaly in transfer ellipse at arrival 
time vs launch date 

8—11(11. Jupiter 70— 71 : Perihelion of transfer orbit vs launch date, 

Type I 

8-11(11). Jupiter 70— 71 : Perihelion of transfer orbit vs launch date, 

Type II 

8—12(1). Jupiter 70— 7 1 : Aphelion of transfer orbit vs launch date, 

Type I 

8—12(11). Jupiter 70— 71 : Aphelion of transfer orbit vs launch date. 

Type II 

8-13(1). Jupiter 70-71 : Inclination of the heliocentric transfer plane 
vs launch date. Type I 

8—13(11). Jupiter 70— 71 : Inclination of the heliocentric transfer plane 
vs launch date. Type II 
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VIII. JUPITER 1970-71: TRAJECTORY PARAMETER GRAPHS (Cont’d) 


Figure 

8-14. 

8-15(1). 

8-15(11). 

8-16. 

8-17. 

8-18(1). 

8-18(11). 

8-19(1). 

8-19(11). 

8 - 20 ( 1 ). 

8 - 20 ( 11 ). 

8 - 21 ( 1 ). 

8 - 21 ( 11 ). 


Jupiter 70— 71: Celestial latitude at arrival time vs launch date 

Jupiter 70— 71 : Asymptotic speed with respect to Jupiter vs 
launch date, Type I 

Jupiter 70— 71 : Asymptotic speed with respect to Jupiter vs 
launch date, Type II 

Jupiter 70— 71 : Angle between incoming zeocentric asymptote 
and arrival planet's orbital plane vs launch date 

Jupiter 70-71: Angle between Jupiter— Sun vector and incoming 
zeocentric asymptote vs launch date 

Jupiter 70— 71 : Declination of zeocentric asymptote vs 
launch date, Type I 

Jupiter 70— 71: Declination of zeocentric asymptote vs 
launch date, Type II 

Jupiter 70-71 : Right ascension of zeocentric asymptote vs 
launch date, Type I 

Jupiter 70-71 : Right ascension of zeocentric asymptote vs 
launch date. Type II 

Jupiter 70— 71 : Angle between planet— Earth vector and 
incoming asymptote vs launch date, Type I 

Jupiter 70— 71 : Angle between planet— Earth vector and 
incoming asymptote vs launch date. Type II 

Jupiter 70— 71: Angle between planet-Canopus vector and 
incoming asymptote vs launch date, Type I 

Jupiter 70— 71: Angle between planet-Canopus vector and 
incoming asymptote vs launch date. Type II 
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Fig. 8—4(1). Jupiter 70— 71; Earth— Jupiter communication distance vs launch date, Type I 
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Fig. 8—4(11). Jupiter 70— 71 : Earth— Jupiter communication distance vs launch date, Type II 














Fig. 8—5(1). Jupiter 70— 71 : Declination of the geocentric asymptote vs launch date, Type I 




































JPL TECHNICAL REPORT NO. 32-77 



190 


DECEMBER JANUARY FEBRUARY MARCH 

LAUNCH DATE 


= 0.95 
= 0.90 


JPL TECHNICAL REPORT NO. 32-77 



Fig. 8—6. Jupiter 70— 71: Right ascension of the geocentric asymptote vs launch date 
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Fig. 8-10. Jupiter 70— 71 : True anomaly in transfer ellipse at arrival time vs launch date 
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Fig. 8—13(11). Jupiter 70— 71 ; Inclination of the heliocentric transfer plane vs launch date, Type II 
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14. Jupiter 70— 71: Celestial latitude at arrival time vs launch date 
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1 5(111. Jupiter 70— 71 : Asymptotic speed with respect to Jupiter vs launch date, Type II 
























16. Jupiter 70-71; Angle between incoming zeocentric asymptote and arrival planet's orbital plane vs launch date 
















1 7. Jupiter 70-71 : Angle between Jupiter— Sun vector and incoming zeocentric asymptote vs launch date 
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Fig. 8—18(1). Jupiter 70— 71: Declination of zeocentric asymptote vs launch date, Type I 
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18(11). Jupiter 70-71 : Declination of zeocentric asymptote vs launch date, Type II 
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t-19(l). Jupiter 70-71 : Right ascension of zeocentric asymptote vs launch date, Type I 
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19(11). Jupiter 70-71 : Right ascension of zeocentric asymptote vs launch date, Type II 
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Fig. 8—20(1). Jupiter 70— 71 : Angle between planet— Earth vector and incoming asymptote vs launch date, Type I 
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Fig. 8— 21 (I). Jupiter 70— 71 : Angle between planet-Canopus vector and incoming asymptote vs launch date, Type I 
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1-21(11). Jupiter 70-71: Angle between planet-Canopus vector and incoming asymptote vs launch date, Type II 
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IX. JUPITER 1972: TRAJECTORY PARAMETER GRAPHS 

figure 

9-1. Jupiter 1972: Minimum injection energy vs launch date 

9-2. Jupiter 1972: Time of flight vs launch date 

9-3. Jupiter 1972: Heliocentric central angle vs launch date 

9—4(1). Jupiter 1972: Earth-Jupiter communication distance vs 
launch date, Type I 

9-4(11). Jupiter 1972: Earth-Jupiter communication distance vs 
launch date, Type II 

9-5. Jupiter 1972: Declination of the geocentric asymptote vs 
launch date 

9—6. Jupiter 1972: Right ascension of geocentric asymptote vs 
launch date 

9—7. Jupiter 1972: Angle betv/een outgoing geocentric asymptote 
and launch planet's orbital plane vs launch date 

9-8(1). Jupiter 1972: Angle between Sun-Earth vector and outgoing 
geocentric asymptote vs launch date, Type I 

9-8(11). Jupiter 1972: Angle between Sun— Earth vector and outgoing 
geocentric asymptote vs launch date. Type II 

9—9(1). Jupiter 1972: True anomaly in transfer ellipse at launch time 
vs launch date, Type I 

9—9(11). Jupiter 1972: True anomaly in transfer ellipse at launch time 
vs launch date. Type II 

9—10. Jupiter 1972: True anomaly in transfer ellipse at arrival time 
vs launch date 

9—11(1). Jupiter 1972: Perihelion of transfer orbit vs launch date, 

Type I 

9—11(11). Jupiter 1972; Perihelion of transfer orbit vs launch date, 

Type II 

9-12(1). Jupiter 1972: Aphelion of transfer orbit vs launch date, Type I 

9-1 2(11). Jupiter 1 972: Aphelion of transfer orbit vs launch date, Type II 

9_13(|), Jupiter 1972: Inclination of the heliocentric transfer plane vs 
launch date. Type I 

9-13(11). Jupiter 1972: Inclination of the heliocentric transfer plane vs 
launch date. Type II 

9-14(1). Jupiter 1972: Celestial latitude at arrival time vs launch date, 
Type I 

9-14(11). Jupiter 1972: Celestial latitude at arrival time vs launch date. 
Type 11 
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IX. JUPITER 1972: TRAJECTORY PARAMETER GRAPHS (Cont’d) 


Figure 

9-1511). 

9-15(11). 

9-16. 

9-17. 

9-18. 

9-19(1). 

9-19(11). 

9-20(1). 

9-20(11). 

9-21(1). 

9-21(11). 


Jupiter 1972: Asymptotic speed with respect to Jupiter vs 
launch date^ Type I 

Jupiter 1972: Asymptotic speed with respect to Jupiter vs 
launch date. Type II 

Jupiter 1972: Angle between incoming zeocentric asymptote 
and arrival planet’s orbital plane vs launch date 

Jupiter 1972: Angle between Jupiter— Sun vector and incoming 
zeocentric asymptote vs launch date 

Jupiter 1972: Declination of zeocentric asymptote vs 
launch date 

Jupiter 1972; Right ascension of zeocentric asymptote vs 
launch date, Type I 

Jupiter 1972: Right ascension of zeocentric asymptote vs 
launch date, Type II 

Jupiter 1972: Angle between planet— Earth vector and incoming 
asymptote vs launch date, Type I 

Jupiter 1972: Angle between planet— Earth vector and incoming 
asymptote vs launch date, Type 11 

Jupiter 1972: Angle between planet— Canopus vector and 
incoming asymptote vs launch date, Type I 

Jupiter 1972: Angie between planet-Canopus vector and 
incoming asymptote vs launch date, Type II 
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Fig. 9“4(ll. Jupiter 1972: Earth— Jupiter communication distance vs launch date, Type I 
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Fig. 9—4(11). Jupiter 1972: Earth— Jupiter communication distance vs launch date, Type II 
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. 9-7. Jupiter 1972: Angle between outgoing geocentric asymptote and launch planet’s orbital plane vs launch date 



















Fig. 9—8(1). Jupiter 1972: Angle between Sun— Earth vector and outgoing geocentric asymptote vs launch date, Type I 
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Fig. 9-8(ll). Jupiter 1972: Angle between Sun-Earth vector and outgoing geocentric asymptote vs launch date, Type II 
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1 1(1). Jupiter 1972: Perihelion of transfer orbit vs launch date, Type I 











1 1(11). Jupiter 1972: Perihelion of transfer orbit vs launch date, Type II 
































Fig. 9—12(1). Jupiter 1972: Aphelion of transfer orbit vs launch date, Type I 























Fig. 9—12(11). Jupiter 1972: Aphelion of transfer orbit vs launch date, Type II 




























13(1). Jupifer 1972: Inclination of the heliocentric transfer plane vs launch date, Type I 
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>—14(11). Jupiter 1972; Celestial latitude at arrival time vs launch date, Type II 
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15(l}. Jupiter 1972: Asymptotic speed with respect to Jupiter vs launch date, Type I 





















16. Jupiter 1972: Angle between incoming zeocentric asymptote and arrival planet’s orbital plane vs launch date 
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17. Jupiter 1972: Angle between Jupiter— Sun vector and incoming zeocentric asymptote vs launch date 
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Fig. 9-20(1). Jupiter 1972: Angle between planet-Earth vector and incoming asymptote vs launch date, Type I 
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Fig. 9—20(11). Jupiter 1972: Angle between planet— Earth vector and incoming asymptote vs launch date, Type II 
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Fig. 9-21(1). Jupiter 1972: Angle between planet-Canopus vector and incoming asymptote vs launch date, Type I 




















'“21(11). Jupiter 1972: Angie between planet-Canopus vector and incoming asymptote vs launch date, Type II 
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X. JUPITER 1973: TRAJECTORY PARAMETER GRAPHS 

figure 

10-1. Jupiter 1973: Minimum injection energy vs launch date 

10-2(1). Jupiter 1973: Time of flight vs launch date, Type I 

10-2(11). Jupiter 1973: Time of flight vs launch date, Type II 

10-3. Jupiter 1973: Heliocentric central angle vs launch date 

10-4(1). Jupiter 1973: Earth-Jupiter communication distance vs 
launch date, Type 1 

10—4(11). Jupiter 1973: Earth-Jupiter communication distance vs 
launch dote, Type II 

10—5(1). Jupiter 1973: Declination of geocentric asymptote vs 
launch date, Type I 

10-5(11). Jupiter 1973: Declination of geocentric asymptote vs 
launch date, Type II 

10-6(1). Jupiter 1973: Right ascension of geocentric asymptote vs 
launch date, Type I 

10—6(11). Jupiter 1973: Right ascension of geocentric asymptote vs 
launch date, Type It 

10-7. Jupiter 1973: Angle between outgoing geocentric asymptote 
and launch planet’s orbital plane vs launch date 

10-8(1). Jupiter 1973; Angle between Sun-Earth vector and outgoing 
geocentric asymptote vs launch date, Type I 

10-8(11). Jupiter 1973; Angle between Sun-Earth vector and outgoing 
geocentric asymptote vs launch date, Type II 

10-9(1). Jupiter 1973: True anomaly in transfer ellipse at launch time 
vs launch date. Type I 

10—9(11). Jupiter 1973: True anomaly In transfer ellipse at launch time 
vs launch date, Type II 

10-10. Jupiter 1973: True anomaly in transfer ellipse at arrival time 
vs launch date 

10—11(1). Jupiter 1973: Perihelion of transfer orbit vs launch date, 

Type I 

10—11(11). Jupiter 1973: Perihelion of transfer orbit vs launch date, 

Type II 

10-12(1). Juplterl973: Aphelion of transfer orbit vs launch date. 

Type I 

10-12(11). Jupiter 1973: Aphelion of transfer orbit vs launch date. 

Type II 
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X, JUPITER 1973: TRAJECTORY PARAMETER GRAPHS (Cont’d) 


Figure 

10-13(1). 

10-13111). 

10-14(1), 

10-14(11). 

10 - 15 ( 1 ). 

10-15(11). 

10-16. 

10-17. 

10-18. 

10-19(11. 

10-19(11). 

10 - 20 ( 1 ). 

10 - 20 ( 11 ). 

10 - 21 ( 1 ). 

10 - 21 ( 11 ). 


Jupiter 1973: Inclination of the heliocentric transfer plane vs 
launch date, Type I 

Jupiter 1973: Inclination of the heliocentric transfer plane vs 
launch date, Type II 

Jupiter 1973: Celestial latitude at arrival time vs launch date, 
Type I 

Jupiter 1973: Celestial latitude at arrival time vs launch date, 
Type II 

Jupiter 1973; Asymptotic speed with respect to Jupiter vs 
launch date, Type I 

Jupiter 1973: Asymptotic speed with respect to Jupiter vs 
launch date. Type II 

Jupiter 1973: Angle between incoming zeocentric asymptote 
and arrival planet’s orbital plane vs launch date 

Jupiter 1973: Angle between Jupiter— Sun vector and incoming 
zeocentric asymptote vs launch date 

Jupiter 1973: Declination of zeocentric asymptote vs 
launch date 

Jupiter 1973: Right ascension of zeocentric asymptote vs 
launch date. Type I 

Jupiter 1973: Right ascension of zeocentric asymptote vs 
launch date, Type II 

Jupiter 1973: Angle between planet— Earth vector and incoming 
asymptote vs launch date, Type I 

Jupiter 1973: Angle between planet— Earth vector and incoming 
asymptote vs launch date, Type II 

Jupiter 1973: Angle between planet-Canopus vector and 
incoming asymptote vs launch date, Type I 

Jupiter 1973: Angle between planet-Canopus vector and 
incoming asymptote vs launch date, Type II 


252 


I 




















Fig. 10—2(1). Jupiter 1973: Time of flight vs launch date, Type I 
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Fig. 10—4(11). Jupiter 1973: Earth— Jupiter communication distance vs launch date, Type II 















Fig. 10-5(1). Jupiter 1973: Declination of geocentric asymptote vs launch date, Type I 




























Fig. 10—5(11). Jupiter 1973: Declination of geocentric asymptote vs launch date, Type II 



























Fig. 10—6(1). Jupiter 1973: Right ascension of geocentric asymptote vs launch date, Type I 





































Fig. 10—6(11). Jupiter 1973: Right ascension of geocentric asymptote vs launch date, Type II 



































Fig. 10—7. Jupiter 1973; Angle between outgoing geocentric asymptote and launch planet's orbital plane vs launch date 









































Pig. 10—8(11). Jupiter 1973: Angle between Sun-Earth vector and outgoing geocentric asymptote vs launch date, Type II 




































Fig. 10—90). Jupiter 1973: True anomaly in transfer ellipse at launch time vs launch date, Type I 





Fig. 10—9(11). Jupiter 1973: True anomaly in transfer ellipse at launch time vs launch date, Type II 
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Fig. 10—12(1). Jupiter 1973: Aphelion of transfer orbit vs launch date, Type I 
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Fig. 10—12(11). Jupiter 1973: Aphelion of transfer orbit vs launch date, Type II 
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Fig. 10—13(1). Jupiter 1973: Inclination of the heliocentric transfer plane vs launch date, Type I 
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LAUNCH DATE 

Fig. 10-13(11). Jupiter 1973: Inclination of the heliocentric transfer plane vs launch date, Type II 
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Fig. 10—14(11). Jupiter 1973: Celestial latitude at arrival time vs launch date, Type II 















Fig. 10—15(1). Jupiter 1973: Asymptotic speed with respect to Jupiter vs launch date, Type I 





















Fig. 10—15(11). Jupiter 1973: Asymptotic speed with respect to Jupiter vs launch date^ Type II 
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Fig. 10—16. Jupiter 1973: Angle between incoming zeocentric asymptote and arrival planet's orbital plane vs launch date 

























JPL TECHNICAL REPORT NO, 32-77 



Fig 10-17. Jupiter 1973: Angle between Jupiter-Sun vector and incoming zeocentric asymptote vs launch date 


















. 10—18. Jupiter 1973: Decimation of zeocentric asymptote vs launch date 
































Fig. 10“19(l). Jupiter 1973: Right ascension of zeocentric asymptote vs launch date, Type I 



















































Fig. 10-19(11). Jupiter 1973: Right ascension of zeocentric asymptote vs launch date, Type II 
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Fig. 10—20(1). Jupiter 1973: Angle between planet— Earth vector and incoming asymptote vs launch date, Type I 































Fig. 10—20(11). Jupiter 1973: Angle between planet— Earth vector and incoming asymptote vs launch date, Type II 
































Fig. 10—21(1). Jupiter 1973: Angle between planet-Canopus vector and incoming asymptote vs launch date, Type I 
























Fig, 10—21(11). Jupiter 1973: Angle between planet-Canopus vector and incoming asymptote vs launch date, Type II 
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NOMENCLATURE 


C , vis viva energy (injection energy of escape 
hyperbola) 

i inclination of heliocentric transfer orbit, deg 

Ra aphelion of transfer orbit, millions of km 

R( Earth-Mercury (or -Jupiter) communication 
distance, millions of km 

Rft perihelion of transfer orbit, millions of km 

T time of flight, days 

V[^ true anomaly in transfer ellipse at launch time, 
deg 

Vp true anomaly in transfer ellipse at arrival time, 
deg 

Vfip asymptotic speed (or hyperbolic-excess speed) 
with respect to Mercury (or Jupiter), km/sec 

celestial latitude of Mercury (or Jupiter) at 
arrival time, deg 

yi^ angle between outgoing geocentric asymptote 
and launch planet’s orbital plane, deg (also 
defined as celestial latitude of outgoing 
asymptote) 

yp angle between incoming hermiocentric^' or 
zeocentric^' asymptote and arrival planet s 
orbital plane, deg 

angle between Sun-Earth vector and outgoing 
geocentric asymptote, deg (also d(^fined as 
Earth-probe-Sun angle) 

(p angle between Mercury-Sun (or Jupiter-Sun) 
vector and incoming hermiocentric (or 
zeocentric) asymptote, deg (also defined as 
Sun-probe-target angle) 

0;, right ascension of hermiocentric (or zeocentric) 
asymptote, deg 

right ascension of geocentric asymptote, deg 

declination of hermiocentric (or zeocentric) 
asymptote, deg 

0,s declination of geocentric asymptote, deg 

^ heliocentric central angle, deg 


^This table of Nomenclature presents only the notation used in tlie 
graphs of Sections IV to X. All other terms are defined at point 
of first mention in Sections I to III. 

^^Hermiocentric refers to Mercury. 

‘ Zeocentric refers to Jupiter. 
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